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ABSTRACT 

The  research  reported  in  this  Final  Technical  Report  may  be  grouped  into  three  main 

areas  as  follows: 

A.  Surface  Guided  Wave  Contributions  to  the  Scattering  of  Sound  by 
Elastic  Shells  in  Water  «  Quantitative  ray  representations  for  the  scattering  by 
shells  should  be  generalizable  to  shells  having  various  shapes  and  give  insight  into 
the  coupling  of  sound  with  the  structure.  The  research  was  concerned  with 
frequencies  above  the  breathing  resonance  or  ring  frequency  of  spheres  and  cylinders 
and  also  with  the  impulse  and  tone  burst  responses.  A  midfrequency  backscattering 
enhancement  near  the  coincidence  frequency  was  observed,  modeled,  and  shown  to 
be  due  to  a  subsonic  wave.  Aspects  of  a  high-frequency  enhancement  were 
investigated  along  with  calculations  of  time-frequency  spectrograms  of  the  impulse 
response, 

B .  Radiation  Mechanisms  for  Evanescent  Tone  Bursts  on  Fluid  Loaded 
Membranes  and  Plates  -  Ordinarily  flexural  waves  with  subsonic  phase 
velocities  on  unbounded  plates  produce  evanescent  wavefields  in  the  surrounding 
fluid  that  are  trapped  to  the  plate.  Inhomogeneities  as  well  as  curvature  cause  energy 
to  be  radiated.  Two  such  mechanisms  were  investigated  experimentally:  (a) 
diffraction  by  a  rigid  barrier  in  the  evanescent  field  of  a  membrane  in  air  and  (b) 
transition  radiation  produced  when  a  wave  packet  travels  down  a  vertical  plate 
through  a  free  water  surface.  Low-frequency  EMAT  transducers  were  used  to 
generate  the  tone  bursts. 

C .  Acoustical  Four-Wave  Mixing  in  Suspensions  -  A  novel  mechanism  was 
investigated  for  the  interaction  of  sound  with  sound  that  relies  on  the  slow  response 
of  the  acoustic  medium  to  the  radiation  pressure  of  a  pump  wave.  Counter 
propagating  800  kHz  pump  waves  establish  a  standing  wave  in  a  saline  suspension  of 
Latex  microspheres.  In  response,  the  spheres  arc  attracted  to  pressure  nodes  and 
form  bands.  The  banded  suspension  acts  as  a  diffraction  grating  for  a  higher 
frequency  probe-wave  tone  burst  giving  rise  to  coherent  scattering  from  the 
suspension  when  the  Bragg  condition  is  satisfied.  Theories  for  the  scattering 
amplitude  using  different  methods  of  analysis  were  found  to  be  mutually  consistent 
and  in  general  agreement  vrith  the  measured  amplitudes  of  steady-state  Bragg 
scattering.  The  temporal  development  of  the  scattering  as  the  grating  initially  forms 
was  also  calculated  and  measured. 
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I .  INTRODUCTION  AND  ORGANIZATION  OF  THIS  REPORT 


This  Final  Technical  Report  is  organized  as  follows.  Section  0  lists  the  personnel 
involved  in  the  research.  Section  in  gives  a  bibliography  of  publications,  reports, 
dissertations,  and  communications  resulting  from  the  sponsored  research.  The  items  are 
listed  to  facilitate  easy  reference  in  the  main  body  of  the  report.  Section  IV  reviews  selected 
research  results  grouped  according  to  the  three  principal  task  areas  originally  proposed;  (A) 
surface  guided  wave  contributions  to  scattering  of  sound  by  elastic  shells  in  water,  (B) 
radiation  mechanisms  for  evanescent  tone  bursts  on  fluid  loaded  membranes  and  plates, 
and  (C)  acoustical  four-wave  mixing  in  suspensions.  Summarized  in  the  final  subsecrion  is 
(D)  supplemental  scattering  research.  Section  V  summarizes  facilities  improvements  that 
facilitated  the  sponsored  research  and  other  ongoing  research.  Section  VI  lists  references 
not  included  in  Section  HI  that  pertain  to  the  discussion  of  resul  ts. 

II.  PERSONNEL 

The  following  persons  participated  in  the  research. 

Graduate  students 

1 .  H.  J.  Simpson:  completed  Ph.D.  in  1992  on  ultrasonic  four-wave  mixing  mediated 
by  a  suspension. 

2.  T.  J.  Matula:  awarded  M.  S.  degree  in  1990  with  project  supported  research. 
Continuing  as  a  Ph.D.  candidate  on  research  pertaining  to  the  radiation  and  .scattering 
by  flexural  waves  on  plates. 

3.  G.  Kaduchak:  awarded  M.S.  degree  in  1992  with  project  related  research.  Working 
as  a  Ph.D.  candidate  (since  1991)  on  scattering  by  elastic  shells. 

4.  D.  H.  Hughes:  Ph.D.  candidate  working  (since  1990)  on  theoretical  and 
computational  problems  pertaining  to  scattering  by  elastic  shells. 

5.  L.  G.  Zhang;  awarded  M.S.  degree  for  project  related  thesis  on  scattering  in  1991. 


6.  S.  P.  Parry:  contributed  to  testing  and  initial  experimentation  in  the  water-tank 
scattering  facility  prior  to  leaving  the  program  for  other  scientific  employment. 

Other  personnel 

1 .  P.  L.  Marston:  principal  investigator 

2.  N.  H.  Sun:  postdoctoral  researcher  who  contributed  to  scattering  theory  and 
computations  in  1990-1991. 

3 .  J.  R.  Filler:  postdoctoral  engineering  support  pertaining  to  the  development  of  the 
scattering  facility  in  1990.  Wliile  on  appointment.  Filler  also  completed  a  publication 
on  previous  research  pertaining  to  fluid-structure  interactions. 

III.  BIBLIOGRAPHY  OF  PUBLICATIONS  AND  REPORTS  FOR  GRANT  N00014-89- 
J-3088 

Research  publications,  pending  publications,  dissertations,  technical  and  internal 
reports,  and  abstracted  oral  presentations  are  listed  below  according  to  category.  Also 
listed  for  reference  in  the  body  of  the  report  is  a  letter  code.  Thus,  for  example,  J 1  refers  to 
the  first  journal  publication.  Each  list  is  roughly  chronological.  Items  supported  in  part  by 
other  sources  are  so  indicated. 

Code  J:  Refereed  Journal  Publications.  Items  that  have  been  accepted  for 
publication  are  so  indicated  while  items  submitted  and  under  review  are  listed  only  as 
"submitted." 

J1 .  S.  G.  Kargl  and  P.  L.  Marston,  "Ray  synthesis  of  the  form  function  for 

backscattering  from  an  elastic  spherical  shell:  Leaky  Lamb  waves  and  longitudinal 
resonances,"  J.  Acoust.  Soc.  Am.  89,  2545-2558  (1991).  [Supported  in  part  by 
N00014-85-C-0141.] 

J2.  J.  R.  Filler,  P.  L.  Marston,  and  W.  C.  Mih,  "Response  of  the  shear  layers  separating 
from  a  circular  cylinder  to  small-amplitude  rotational  oscillations,"  Journal  of  Fluid 
Mechanics  231, 481-499  (1991).  [Supported  in  part  by  N()(X)14-85-C-0141.] 


J3.  N.  H.  Sun  and  P.  L.  Marston,  "Ray  synthesis  of  leaky  Lamb  wave  contributions  to 
backscattering  from  thick  cylindrical  shells,"  J.  Acoust.  Soc.  Am.  91, 1398-1402 
(1992). 

J4.  L.  G.  Zhang,  N.  H.  Sun,  and  P.  L.  Marston,  "Midfrequency  enhancement  of  the 
backscattering  of  tone  bursts  by  thin  spherical  shells,"  J.  Acoust.  Soc.  Am.  91, 
1862-1874(1992). 

J5.  C.  K.  Frederickson  and  P.  L.  Marston,  "Transverse  cusp  diffraction  catastrophes 
produced  by  the  reflection  of  ultrasonic  tone  bursts  from  a  curved  surface  in  water: 
Observations,"  J.  Acoust.  Soc.  Am.  92,  2869-2877  (1992).  [Supported  in  part  by 
N00014-85-C-0141.] 

J6.  P.  L.  Marston  and  N.  H.  Sun,  "Resonance  and  interference  scattering  near  the 
coincidence  frequency  of  a  thin  spherical  shell:  an  approximate  ray  synthesis."  J. 
Acoust.  Soc.  Am.  (accepted  for  publication). 

J7 .  G.  Kaduchak  and  P.  L.  Marston,  "Observation  of  the  midfrequency  enhancement  of 
tone  bursts  backscattered  by  a  thin  spherical  shell  in  water  near  the  coincidence 
frequency,"  J.  Acoust.  Soc.  Am.  (accepted  for  publication). 

J8.  T.  J.  Matula  and  P.  L.  Marston,  "Electromagnetic  acoustic  wave  transducer  for  the 
generation  of  acoustic  evanescent  waves  on  membranes  and  an  optica!  wave-number 
selective  detector,"  J.  Acoust.  Soc.  Am,  (submitted  in  1992). 

J9.  T.  J.  Matula  and  P.  L.  Marston,  Diffraction  of  evanescent  wave  tone  bursts  on  a 
membrane  in  air,"  J.  Acoust  Soc.  Am.  (accepted  for  publication). 

Code  B:  Chapters  in  Books. 

B 1 .  P.  L.  Marston,  S.  G.  Kargl,  and  N.  H.  Sun,  "Elastic  resonance  amplitudes  described 
by  generalized  GTD  and  by  product  expansions  of  the  S -matrix,"  in  Acoustic 
Resonance  Scattering  edited  by  H.  Uberall  (Gordon  and  Breach,  New  York,  1992) 
pp.  305-333.  [Supported  in  part  by  N0(X)14-85-C-0141.] 
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B2.  P.  L.  Marston,  "Geometrical  and  Catastrophe  Optics  Methods  in  Scattering,"  in  High 
Frequency  and  Pulsed  Scattering.  Physical  Acoustics,  Vol.  XXI,  edited  by  R.  N. 
Thurston  and  A.  D.  Pierce  (Academic  Press,  Boston,  1992)  pp.  1-234.  [Supported 
in  part  by  N00014-85-C-0141,] 

B3.  P.  L.  Marston,  "Quantitative  Ray  Methods  for  Scattering,"  Handbook  of  Acoustics, 
edited  by  M.  J.  Crocker,  (John  Wiley  Press),  submitted. 

B4.  P.  L.  Marston,  "Introductory  Chapter  -  Ultrasonics,  Quantum  Acoustics  and  Physical 
Effects  of  Sound,"  Handbook  of  Acoustics,  edited  by  M.  J.  Crocker,  (John  Wiley 
Press),  submitted. 

Code  T:  Thesis,  Dissertations,  or  Technical  Reports. 

T1 .  Ligang  Zhang,  "Scattering  of  tone  burst  from  spherical  shells:  computations  based 
on  Fourier  transform  method,"  M.S.  thesis.  Department  of  Physics,  Washington 
State  University  (1991),  154  pages.  Grant  Technical  Report  Number  1;  DTIC  report 
number  A240642. 

T2.  Harry  J.  Simpson,  "Interaction  of  Sound  with  Sound  by  Novel  Mechanisms: 

Ultrasonic  four- Wave  Mixing  Mediated  by  a  suspension  and  Ultrasonic  Three-Wave 
Mixing  at  a  Free  Surface,"  Ph.D.  dissertation.  Department  of  Physics,  Washington 
State  University  (1992),  318  pages.  Grant  Technical  Report  Number  2;  DTIC  report 
number  A254657. 

Code  P:  Papers  In  Conference  Proceedings. 

PI.  P.  L.  Marston  and  S.  G.  Kargl,  "Scattering  from  hollow  shells:  Quantitative  ray 

representations  of  amplitudes,"  in  Proceedings  of  the  International  Congress  on  Recent 
Developments  in  Air  and  Structure  Bourne  Sound  and  Vibration,  edited  by  M.  J. 
Crocker  (Auburn  University,  1990)  pp.  565-69.  [Supported  in  part  by  N00014-85-C- 
0141.] 

P2.  Invited  paper:  P.  L.  Marston,  L.  Zhang,  N.  Sun,  G.  Kaduchak,  D.  H.  Hughes,  "Ray 
representations  of  the  backscattering  of  tone  bursts  by  shells  in  water:  Calculations 
and  related  experiments,"  in  Proceedings  of  the  Second  International  congress  on 


Recent  Developments  in  Air  -  and  Structure-Bome  Sound  and  Vibration,  edited  by  M. 
J.  Crocker  and  P.  K.  Raju  (Auburn  University,  1992),  pp.  1203-1210. 

P3.  L.  Zhang,  N.  Sun,  and  P.  L.  Marston,  "Mid-frequency  enhancement  of  the 

backscattering  of  tone  bursts  by  thin  spherical  shells:  ray  approximation  and  exact 
Fourier  synthesis,"  in  Proceedings  of  the  14th  International  Congress  on  Acoustics, 
edited  by  Li  Peizi  (Beijing,  1992),  pp.  A-7-9.1, 2. 

Code  I;  Internal  Reports.  These  include  various  reports  on  file  that  were  not  widely 

distributed.  Copies  will  be  available  for  a  limited  time  form  the  principal  investigator. 

Some  of  these  have  been  superseded  by  other  manuscripts  that  were  submitted  for 

publication. 

11.  T.  J.  Matula,  "Electromagnetic  Generation  of  Plane  Waves  on  Membranes  for  the 
Production  of  Acoustic  Evanescent  Waves,"  M.S.  degree  project  report.  Department 
of  Physics,  Washington  State  University  (1990). 

12.  T.  J.  Matula,  "Diffraction  of  evanescent  wave  tone  bursts  on  a  membrane  in  air  and 
related  observations  of  bursts  guided  by  plates  in  water,"  March,  1992,  Note:  This 
paper  (and  associated  presentation)  by  Matula  resulted  in  the  Second  Place  Award  for  a 
student  paper  on  Structural  Acoustics  and  Vibration,  Salt  Lake  City  Meeting  of  the 
Acoustical  Society  of  America,  April,  1992. 

13.  G.  Kaduchak,  "Observation  of  the  midfrequency  enhancement  of  tone  bursts 
backscattered  by  a  thin  spherical  shell  in  water,"  March,  1992. 

14.  D.  H,  Hughes,  "Impulse  response  for  backscattering  by  a  thin  spherical  shell 
displayed  in  time  and  frequency,"  March,  1992. 

15.  G.  Kaduchak,  "Observation  of  the  E^  Diffraction  Catastrophe  by  Optical  Scattering 
from  Levitated  Drops  in  the  Rainbow  Region,"  M.S.  degree  project  report. 

Department  of  Physics,  Washington  State  University.  [Supported  in  part  by  NOOOH- 
85-C-0141.] 

Code  M:  Oral  Presentations  at  Professional  Meetings.  Presentations  at 
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government  research  centers  or  universities  are  not  listed. 


Ml.  P.  L.  Marston,  "Applications  to  radiation  and  scattering  of  a  wavefront  construction 
using  the  involute  of  circular  virtual  caustics,"  J.  Acoust.  Soc.  Am.  87, 50  (1990). 

M2.  S.  G.  Kargl  and  P.  L.  Marston,  "Anomalies  in  the  ray  synthesis  of  backscattering 
from  hollow  elastic  spherical  shells,"  J.  Acoust.  Soc.  Am.  88, 15  (1990). 
[Supported  in  part  by  N00014-85-C-0141.] 
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IV  RESEARCH  PROJECTS 


This  discussion  summarizes  the  research  results  in  the  context  of  previous  work  by 
others  and  by  Marston  and  associates.  The  reader  should  see  the  various  publications  and 
reports  listed  in  Section  m  for  theoretical  and  technical  details. 

A.  Surface  Guided  ^ave  Contributions  to  the  Scattering  of  Sound  by 
Elastic  Shells  in  Water. 

1.  Ray  Methods  for  Leaky  Waves  on  Spheres  and  Cylinders. 

An  important  advance  in  the  development  of  ray  methods  for  scattering  by  shells  was 
the  realization  that  the  coupling  of  the  acoustic  field  with  guided  wave  amplitudes  could  be 
expressed  in  a  simple  form.^  The  application  of  that  analysis  to  leaky  Lamb  waves  on 
spherical  shells  was  confirmed  with  tone  burst  measurements  by  Kargl  and  Marston.^  The 
relevant  coupling  coefficient  G/  for  the  Ith  class  of  guided  waves  was  expressed  in  terms  of 
the  radiation  damping  of  the  guided  wave  p/  and  geometric  factors  associated  with  the  location 
of  the  internal  virtual  caustics  from  which  the  outgoing  wave  appeared  to  radiate.*  The  ray 
method  also  facilitates  the  identification  of  elastic  contributions  by  their  arrival  time  when 
short  tone  bursts  are  used.  Extensions  should  be  possible  for  shells  with  slowly  varying 
curvature  and  internal  structure.  Furthermore,  since  G/  appears  to  be  well  approximated 
directly  in  terms  of  guided  wave  parameters,  it  becomes  unnecessary  to  analytically  and 
numerically  evaluate  it  from  the  Sommerfeld-Watson  transform  as  was  done  by  Williams  and 
Marston^  for  Rayleigh  waves  on  solid  spheres.  Kargl  and  Marston'^'^Jl  completed  other 
computations,  tests,  and  extensions  of  the  method  for  forward  and  backward  scattering  by 
spherical  shells.  While  the  emphasis  was  on  resonances  associated  with  leaky  Lamb  waves, 
the  earlier  experiments^  were  also  extended  to  forward  scatteiing^  and  some  aspects  of 
thickness  resonance  weie  modeled.^***  In  related  work  a  ray  model  for  the  backscattering 
contribution  due  to  repeated  reverberation  (or  bouncing)  of  normally  incident  longitudinal 
waves  was  developed  and  tested  by  comparison  with  exact  partial-wave  series  calculations  for 
a  hollow  fluid  shell.  This  was  first  achieved  by  Kargl  and  Marston^  for  a  sphere  and  was 
later  extended  with  support  from  this  grant  by  Marston  and  Sun  (Section  IV  D)  for  a  circular 
cylinder. 

An  important  result  from  this  grant  concerns  backscattering  conuibutions  by 
subsonic  waves  on  shells.  It  is  appropriate  to  review  the  ray  picture  for  leaky  and  subsonic 
cases.  Figure  1  illustrates  the  usual  ray  picture  for  backscattering  contributions  by  the  /th 


12 


guided  wave  where  the  phase  velocity  of  the  wave  along  the  outer  surface  of  the  shell  c/ 
exceeds  the  speed  of  sound  c  in  the  surrounding  water.  Such  supersonic  waves  usually 
continuously  leak  or  radiate  sound  at  an  angle  9/ =arcsin  (c/c/)  relative  to  the  local  shell 
normal  as  determined  by  the  trace  velocity  matching  condition  (see  however,  discussion 
below).  To  facilitate  the  extension  of  midfrequency  ray  methods  to  thin  shells,  it  is 
necessary  to  allow  for  subsonic  guided  waves  where  c/  <  c.  Figure  2  shows  the  modified 
ray  picture  where  the  coupling  is  through  an  evanescent  region  having  a  thickness  (b/ 
-  a)  =  (uc/c/)(l  -  c//c)  that  diverges  as  the  outer  radius  a  of  the  shell  diverges.  Apparent 
from  truncations  or  other  inhomogeneities,  the  subsonic  waves  considered  would  not 
radiate  in  this  flat  plate  limit  and  are  sometimes  described  as  "trapped".  The  existence  of 
such  waves  for  plates  in  water  has  long  been  known.^'l  ^  The  radiation  and  coupling  with 
acoustic  fields  introduced  by  curvature  is  analogous  with  the  bending  losses  of  fiber-optic 
waveguides.®^  The  sponsored  analysis  and  experiments  (Publications  J4,  J6,  J7,  Tl,  P2, 
and  P3)  give  quantitative  scattering  amplitudes  based  on  the  moditied  ray  picture  in  Fig.  2 
and  are  summarized  below. 

2.  Midfrequency  Enhancement  of  the  Backscattering  of  Tone  Bursts  Near  the  Coincidence 
Frequency. 

The  contribution  of  a  subsonic  guided  wave  to  the  backscattering  by  thin  spherical 
shells  is  especially  significant  when  ka  is  in  the  general  vicinity  of  a/h  where  h  =  fl  -  b  is  the 
shell  thickness  and  k  »  oVc.  [This  corresponds  to  kh  close  to  unity  which  for  a  metallic 
plate  or  shell  in  water  is  close  to  the  "coincidence  condition."  Here  the  phase  velocity  of 
the  flexural  wave  (without  fluid  loading)  coincides  to  that  of  sound  in  the  surrounding 
water.  Some  authors  call  this  the  critical  coincidence  frequency.]  For  the  purpose  of 
illustrating  a  ray  model  for  such  contributions,^^  consider  the  example  of  a  stainless  steel 
304  shell  with  a/h  =  40.  The  relevant  guided  wave  properties  are  computed  by  finding  the 
complex  V  roots  of  £)v(ka)  =  0  as  previously  described2.4,3  where  Dn  is  the  denominator  of 
the  nth  term  of  the  exact  partial  wave  series  for  the  form  function.  (The  method  is  based  on 
the  Watson  transform.)  The  resulting  properties  for  the  case  under  consideration  are 
shown  in  Fig.  3.  The  wave  labeled  with  long  and  short  dashes,  designated  by  /  =  ao-,  is 
subsonic  throughout  the  region  of  interest  but  displays  a  noticeable  rise  in  the  radiation 
damping  parameter  P/  as  c/  approaches  c.  This  rise  in  coupling  appears  to  be  associated 
with  the  decrease  in  the  thickness  of  the  evanescent  region  near  coincidence^^*^^  and  causes 
an  enhanced  backscattering  as  may  be  seen  from  the  ray  analysis  summarized  below. 
Extension  of  the  analysis  of  guided  wave  scattering  contributions  by  distinct  leaky 
waves  to  the  present  case  of  subsonic  waves  is  possible  by  appreciating  the  importance 


0  20  40  M  to  100 


ka 

Figure  3  -  Nomialized  phase  velocity  (a)  and  radiation  dating  (b)  of  surface  ^ded  waves  on  a 
3M  stainless  steel  shell  in  water  with  h/a>  0.025.  The  curves  are  identified  as  follows:  long  and 
short  dashes  (/  ■  ao.).  dashes  (/  •  a(0,  and  solid  (/  ■  ao). 
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Rguit  4 -The  solid  curve  fives  die  ray  model  in  Eq.  (1)  for  the  •o.fnided  wave  edioaiiqilitude 
die  region  of  the  enhancement  neffka*  46.  Thepoibts  are  from  the  eotact  echo  anqilitude 
detomined  by  a  Fourier  synthesis  as  iUuatiated  in  I^.  5. 


of  the  caustic  radius  b/  to  the  backscattering  amplitude.^^  The  resulting  contribution  of  the 
mth  circumnavigation  to  the  form  function  f  has  the  magnitude 

If  I »  StcP  (c/c)  exp(-7tp  -  2jcmp ),  (1) 

m/  /  I  /  / 

where  for  m  =  0,  the  wave  has  traveled  only  around  the  backside  of  the  sphere  as  shown  in 
Fig.  2.  The  solid  curve  in  Fig.  4  gives  lfo/1  for  /  =  ao-  while  the  dashed  curve  omits  the 
factor  c/c/ ;  Ifo/I  for  the  ao-  wave  shows  a  strong  enhancement  that  peaks  near  ka  =  46.3 
where  the  value  of  Ifo/I «  3.17.  It  may  be  shown^^  that  the  value  of  Ifo/I  at  such  a  peak 
depends  only  weakly  on  the  shell  thickness  and  material  parameters  since  the  ka 
dependance  is  dominated  by  that  of  P/.  The  peak  Ifo/I  occurs  close  to  where  P/  is  l/ji. 

Equation  (1)  is  directly  applicable  to  the  calculation  of  the  enhanced  backscattering 
by  tone  bursts  that  are  sufficiently  short  that  the  echoes  associated  with  different  values  of 
m  don't  overlap.  The  incident  burst  must  be  sufficiently  long  that  the  effects  of  dispersion 
are  weak.  The  time  domain  echoes  for  such  tone  bursts  were  calculated  by  a  Fourier 
synthesis  that  used  the  exact  partial  wave  series.'*'*’^*  This  was  done  for  several  bursts  of 
different  ka.  Figure  5  shows  the  resulting  calculated  amplitude  for  a  20  cycle  burst  with  ka 
=  46.  The  dimensionless  time  units  are  T  =  ct/a.  The  normalization  is  such  that  specular 
reflection  by  a  fixed  rigid  sphere  of  radius  a  has  a  unit  amplitude.  The  earliest  contribution 
from  the  shell  in  Fig.  5  is  a  specular  reflection  of  close  to  unit  amplitude.  This  is  followed 
by  the  m  =  0  and  m  =  1  contributions  of  the  /  =  ao-  wave  where  the  amplitude  of  the  m  =  0 
echo  of  3.07  is  close  to  the  aforementioned  predicted  enhancement.  The  identity  of  these 
echoes  was  also  confirmed  by  comparison  of  their  arrival  times  with  ray  theory.  The 
amplitude  of  the  m  =  0,  /  =  ao.  echo  was  similarly  determined  for  several  other  values  of  ka 
and  are  plotted  as  the  points  in  Hg.  4.  The  comparison  with  Eq.  (2)  gives  strong  support 
for  the  ray  model  of  the  enhancement.  It  is  noteworthy  that  the  later  (m  =  1)  echo  in  Fig.  5 
is  weaker  by  a  factor  clo.se  to  the  expected  value  of  exp(-27iP/)  «  exp(-2) »  0.135. 
Furthermore,  the  predicted  amplitudes  are  weak  for  the  distinct  ray  contributions  of  the  /  = 
ao  and  sq  guided  waves  whose  properties  are  also  shown  in  Fig.  3.  Consequently  those 
contributions  only  weakly  affect  the  time  record  in  Fig.  5.  The  enhanced  contribution  to 
the  backscattering  by  the  ao-  wave  may  be  especially  useful  for  inverse  problems. The 
echo  is  generally  much  larger  than  leaky  wave  echoes  previously  studied.2.3 

3.  Laboratory  Observations  of  the  Midfirequency  Enhancement. 

The  existence  and  general  magnitude  of  the  enhancement  described  above  was 
confirmed  with  experiments  in  a  large  redwood  tank  of  water.J^  [See  Section  V  for  a 
description  of  the  facility.]  The  volume  of  water  available,  approximately  6500  gallons. 
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(a)  f=340kHz  ^  (c)  fco=55.0 


Figure  8  ~  (a)  and  (b)  represent  expoimental  time  records  like  in  Fig.  6{a)  except  the  carrier  frequency  is  raised, 
(c)  and  correspond  to  the  respective  Fourier  synthesis  calculations. 


was  such  that  reflections  from  the  sides,  top,  and  bottom  of  the  tank  did  not  occur  in  the 
time  region  of  interest  The  target  was  a  stainless  steel  304  shell  with  a  radius  a  ==  38.1  mm 
and  a  radius-to-thickness  ratio  fl/h=  43.8.  The  predicted  enhancement  is  for  frequencies 
near  310  kHz.  Figure  6(a)  shows  a  representative  time  recoid  in  the  region  of  greatest 
enhancement  of  the  guided  wave  echo  which  is  the  second  of  the  echoes.  This  record, 
which  is  for  320  kHz  manifests  the  important  features  of  the  Fourier  synthesis  for  this  shell 
shown  in  Fig.  6(b).  A  substantial  shift  of  the  carrier  frequency  was  observed  to  cause  a 
reduction  in  the  guided  wave  echo  magnitude  relative  to  that  of  the  specular  reflection  and 
the  observed  behavior  is  in  at  least  qualitative  agreement  with  the  shape  of  Ifo/I  as  is  evident 
from  the  comparison  shown  in  Fig.  7.  The  results  are  supportive  of  the  ray  model.^7 
Figuies  8  and  9  show  how  the  echoes  are  in  detailed  agreement  with  predictions  at  other 
frequencies. 


4.  Midftequency  Enhancement  of  the  Form  Function. 

The  ray  model  was  also  used  to  synthesize  the  exact  form  function  f  for  steady-state 
backscattering  in  the  region  of  the  midfrequency  enhancement  for  the  spherical  shell 
considered  in  Fig.  3  with  a/h  =  40.  The  dashed  curve  in  Fig.  10  shows  the  Ifl  from  the 
exact  partial  wave  series  and  the  solid  curve  shows  Ifrayl  where^^ 


^rav  “  ^SD 

ray  sp  a^  a^ 


(2) 


where  fsp  is  the  specular  contribution  (here  neglecting  a  smaU  curvature  correction)  and  f/ 
for  /  =  SQ  and  ao  are  leaky  wave  contributions  computed  as  discussed  previously. In  the 
region  shown  If/I  is  negligible  for  /  =  ao  but  for  /  s  sq  it  gives  rise  to  the  vary  narrow 
resonance  spikes.  The  enhancement  is  associated  with  the  /  =  ao-  term  which  is  computed 
by  summing  ray  terms  of  the  form  fmi  as  in  Eq.  (1)  but  including  phase  information.  The 
result  reduces  to 


-G^exp[-ttp^->-i-nJ 
1  +  exp(-2jtp^  +  i2jtxc/c^) 


(3,4) 


where  11/  =  wxc/c/  -  (7i/2),  x  =  kn,  and  except  for  the  expression  for  arg(G/)  =  <|>/,  the 
generalization  to  this  case  of  a  subsonic  wave  follows  from  the  ray  geometry  in  Fig.  2  and 
related  considerations.  ^  >14  The  approximation  for  used  in  Fig.  10  was  obtained  by 
assuming  that  (|>/  varies  linearly  with  c/  as  c/  approaches  c. 


(Rc  <c^<c)«7c[(3/2)-(2F/3)],  F  =  [(c/c)- v]/(1  -  V)  ,  (5,6) 

but  that  (j)/  becomes  constant  for  c//c  below  a  particular  value  of  V  <  1  so  that  d^ici  <  Vc)  = 
3n/2.  Note  that  F  varies  from  0  to  1  in  the  domain  of  Eq.  (5).  Results  of  Ho  and  Felsen 
were  used  as  a  guide  for  formulating  the  limiting  values  of  ((»/  though  the  value  of  V  =  0.91 
used  for  the  synthesis  in  Fig.  10  was  determined  empirically.  Inspection  of  Fig.  3(b) 
shows  that  for  ka  near  40,  (3/  is  sufficiently  small  for  the  ao-  wave  that  contributions  from 
repeated  circumnavigations  are  important  giving  rise  to  a  distinct  resonance  structure.  With 
increasing  ka  there  is  an  enhancement  in  If/I  associated  primarily  with  the  m  =  0  contribution 
plotted  in  Fig.  4.  For  ka  ^  50,  P/  is  so  large  that  repeated  circumnavigations  are  no  longer 

important  and  the  gradual  periodic  variations  in  Ifl  shown  in  Fig.  10  are  due  to  the 
interference  with  the  specular  reflection.  This  was  verified  with  a  simple  ray  calculation  of 
the  quasiperiod  Aka  »  2n:[2  +  (itc/c/)]**  *  1.2.  A  ray  synthesis  was  also  demonstrated  for 
backscattering  from  a  circular  cylindrical  stainless  steel  shell  also  with  a/h  =  40  and  the 
agreement  with  the  exact  result  was  slightly  improved  over  that  of  Fig.  10.^"^  The 
magnitude  of  the  enhancement  relative  to  the  strength  of  the  specular  contribution  is  weaker 
for  cylinders  than  for  spheres  due  to  the  absence  of  axial  focusing  in  the  cylinder  case.* 

It  is  noteworthy  that  there  has  been  considerable  attention  by  other  groups  given  to 
the  enhancement  of  Ifl  near  coincidence  as  calculated  directly  from  the  exact  partial-wave 
series  as  in  the  dashed  curve  in  Fig.  10.  (See  e.g.  Refs.  12-14  and  other  publications  cited 
therein.)  The  approach  usually  taken  to  understand  the  broad  enhancement  is  to  first  locate 
resonances  in  the  complex  frequency  plane  for  different  partial  waves  n.  From  the  trajectory 
of  resonances  an  attempt  is  made  to  infer  the  existence  and  phase  velocities  of  the  various 
surface  guided  waves.  It  is  noteworthy  that  phase  velocity  curves  found  from  the  complex 
v  roots  such  as  shown  in  Fig.  3(a)  are  consistent  with  results  for  a  flat  plate  near  the 
coincidence  frequency  and  are  perhaps  clearer  than  some  of  the  published  dispersion  curves 
based  on  complex  frequency  plane  roots.  (The  results  of  Sammelmann  et  al.J2  exhibit  the 
same  behavior  for  c/  as  Fig.  3(a).)  The  point  here  is  that  in  the  region  of  the  coincidence 
anomaly,  the  locus  of  complex-v  roots  used  in  ray  theory  are  easy  to  interpret  since  P/(ka)  is 
given  directly.  Another  advantage  of  complex-v  plane  ray  theory  is  evident  from  the 
simplicity  of  Eqs.  (2)-(4).  The  ray  approach  provides  a  natural  mathematical  basis  for 
describing  the  wavefield  in  the  region  of  rapidly  varying  damping  (3/  (Fig.  3(b))  and  that,  at 
least  in  principle,  the  ray  approach  can  be  generalized  to  other  shapes  for  which  no  partial 
wave  series  exists.  There  are  advantages  of  ray  approaches  based  on  complex-v  roots  both 
for  the  time  domain  envelopes.  Fig.  4,  and  for  the  frequency  responses.  Fig.  10. 


5.  Impulse  Response  of  a  Thin  Spherical  Shell  and  a  Significant  Wave  Packet  at  the 
Coincidence  Frequency. 

It  is  well  known  that  for  linear  systems,  the  response  to  an  arbitrary  "input"  may  be 
expressed  as  a  convolution  of  the  input  with  the  system's  impulse  response.  It  follows  that 
an  understanding  and  analysis  of  a  shell's  farfield  scattering  in  response  to  a  5-function 
pressure  impulse  will  give  insight  into  the  scattered  pressure  for  an  arbitrary  incident  wave. 
This  method  of  approach  to  scattering  from  shells  appears  to  be  underutilized.  Figure  1 1(a) 
shows  the  farfield  impulse  response  for  the  thin  shell  considered  in  Figs.  2-5  computed  by 
Zhang,  Sun,  and  Marston^^  from  a  numerical  Fourier  transform  of  the  exact  partial  wave 
series.  The  size  of  the  initial  specular  spike  and  the  relative  normalization  to  the  subsequent 
guided  wave  contributions  are  discused  in  Ref.  J4.  What  is  important  to  the  present 
discussion  is  the  prominent  feature  labeled  ao.  that  has  a  frequency  close  to  the  coincidence 
frequency  of  the  shell.  The  other  guided  wave  contributions  are  lower  in  amplitude  and 
relatively  chirped.  The  large  magnitude  of  the  ao-  guided  wave  echo  in  Fig.  5  can  be 
understood  as  the  convolution  of  that  wave  packet  with  the  incident  tone  burst.  A  weaker 
subsequent  (m  =;  1)  wave  packet  for  the  ao-  wave  yields  the  lower  amplitude  late  echo  in 
Fig.  5. 

It  is  worth  considering  the  earliest  wave  packet's  large  contribution  to  the  impulse 
response  from  an  approximate  analytical  Fourier  transform  of  the  associated  contribution  to 
the  form  function  as  predicted  by  ray  theory.  This  is  because  such  a  calculation  is  a 
sensitive  test  of  ray  methods  and  may  in  principle  be  extended  to  complicated  objects.  (The 
initial  specular  spike  may  be  predicted  from  ray  theory  as  discussed  by  Zhang,  Sun,  and 
Marston.J^)  To  obtain  an  approximate  analysis  of  the  ao-  wave  packet  Marston^^l2  used  a 
stationary  phase  approximation  to  analytically  approximate  the  Fourier  transform  of  the 
relevant  spectral  contribution,  fT,  =  qj.  First  it  was  noted  that  Ifm  =s  o,/l  as  shown  in  Fig.  4 
could  be  well  approximated  by  a  Gaussian  function  centered  at  the  peak  ka  =  x/ «  46.4. 

The  width  of  the  Gaussian  was  expressed  in  terms  of  ray  parameters  evaluated  at  that  ka. 
Then  the  approximate  phase  from  Eqs.  (5)  and  (6)  and  the  propagation  delay  were 
introduced.  Only  the  leading  effect  of  dispersion  (the  dependence  of  c/  on  ka)  could  be 
introduced  and  certain  other  simplifications  were  made  which  led  to  a  narrow  range  of 
possible  parameters  for  the  scattered  wave  packet.  The  resulting  approximate  analytical 
wave  packet  mav  be  reduced  to  the  simple  form  for  the  ao-  wave 

x=x,  +  (AT/T'),  (7) 

where  the  time  shift  relative  to  the  center  of  the  packet  at  T/  is  AT  =  T  -  T/,  The  amplitude 
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Figure  12  ••  Time*  frequency  analysis  of  the  impulse  response  for  backscattering  by  the 
2.5%  thick  shell  considered  in  Figs.  3-5.  The  dynamic  range  is  such  that  the  digest 
regions  correspond  to  signals  that  are  orders-of-magnitude  larger  than  the  light  regions. 

The  earliest  and  strongest  of  a  sequence  of  ao-  wave  packets  is  labeled  as  is  the  prompt 
large  thickness-resonance/si-wave-threshold  feature,  both  of  which  are  stronger  in  this  plot 
than  the  specular  reflection.  Features  associated  with  circumnavigations  of  dispersive  leaky 
Lamb  waves  are  also  evident. 
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Figure  13  *•  Plot  of  the  Gaussian  smoothed  Wigner  distribution  function  (SWDF)  of  the 
impulse  response  for  the  same  shell  as  Hg.  12  wt  for  x  ■  ka  from  0  to  100.  The  SWDF  is 
normalized  to  the  maximum  value  which  is  the  peak  from  the  ao.  wave  packet.  The  ridge  at 
T  «  1  is  due  to  the  specular  reflection. 


factor  A  is  predicted  by  the  analysis  as  are  estimates  for  T/,  T',  and  'P.  Figure  11(b) 
compares  an  exact  impulse  response  (the  solid  curve)  on  an  expanded  scale  with  the 
approximate  analytical  form  given  by  Eq.  (7).  The  exact  response  was  computed  as  in  Fig. 
1 1(a)  but  with  a  different  offset  Tq  so  that  the  dimensionless  time  T  is  relative  to  the 
specular  spike  that  is  outside  of  the  range  of  Fig.  1 1(b).  Also  shown  is  the  envelope  of  the 
wave  packet  as  predicted  by  ray  theory.  The  important  result  is  that  the  simple  ray 
approximation  describes  moderately  well  the  exact  packet.  One  shortcoming  of  the  theory 
(in  addition  to  the  imprecision  of  the  wave  packet  parameters)  is  that  the  exact  packet 
displays  a  weak  low-ffequency  to  high-frequency  chirp  not  analyzed  for  the  degree  of 
dispersion  retained  in  the  approximation.  (The  wave  packet  in  Eq.  (7)  has  a  fixed 
frequency.)  Nevertheless,  there  is  no  doubt  as  to  the  basic  cause  of  this  feature  of  the 
impulse  response. 

6.  Time-Frequency  Analysis  of  Backscattering  by  Spherical  Shells  and  Applications. 

Yen  et  al.l^  explored  the  use  of  a  modified  version  of  the  Wigner  distribution 
function  (WDF)  for  displaying  a  time-frequency  analysis  of  transients  scattered  by  elastic 
objects.  In  our  research  D.  H.  Hughes  applied  methods  of  smoothing  the  WDF  discussed 
by  Nuttall*^  to  scattering  by  spherical  shells.  The  example  considered  below  is  that  of 
backscattering  by  the  same  shell  considered  in  Figs.  3-5.  The  response  to  an  impulse  is 
computed  for  the  time  domain  in  Fig.  1 1(a)  where  it  is  shown  that,  in  addition  to  the 
specular  reflection,  there  is  a  large  wave  packet  associated  with  the  /  =ao.  wave  with  the  m 
=  0  ray  path  shown  in  Fig.  2.  The  impulse  response  in  the  time-frequency  domain 
displays  the  spectral  evolution  of  this  and  other  wave  packets.  Figure  12  shows  the 
smoothed  WDF  for  the  impulse  response  based  on  the  exact  f  from  the  partial-wave 
series.^'^*^’*^  The  axes  are  in  dimensionless  kfl  and  T  =  tc/a  wheie  the  time  offset  used 
is  the  same  as  that  of  Fig.  5  so  that  the  earliest  contri  >ution  is  the  specular  /eflection  at  T  = 

1  in  the  narrow  time  resolution  limit  Various  contributions  can  be  identified  by  comparing 
the  time  of  the  conhibution  with  predictions  from  ray  theory  as  discussed  e.g.  in  Ref.  2  and 
J4.  The  dark  bar  centered  on  T  =  1  is  the  specular  reflection  and  the  dark  patch  centered 
near  T  •>>  S.S  and  ka  »  46  is  from  the  /  =  ao.,  m  =  0  wave  packet  associated  with  the 
midfrequency  enhancement  The  central  frequency  is  close  to  the  coincidence  frequency 
and  the  timing  of  the  packet  follows  from  the  ray  analysis.  Leaky  wave  contributions  that 
can  be  easily  identified  include  those  due  to  the  /  =  so  and  ai  generalizations  of  Lamb  waves 
to  the  present  case  of  a  hollow  fluid-loaded  shell.  Figure  13  shows  another  way  of  plotting 
the  smoothed  WDF  that  also  reveals  the  features  discussed  above. 
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Figure  12  shows  that  there  is  a  large  prompt  contribution  to  the  backscattering  for 
ka  in  the  region  from  approximately  455  to  485.  'fhis  enhancement  is  in  the  vicinity  of  the 
threshold  for  the  si  Lamb  wave  as  well  as  the  longitudinal  resonance  discussed  in  Ref.  6 
and  J1  that,  for  the  present  shell  [Eq.  (9),  below],  is  at  (ka)LR  =  482.  Since  Fig.  12  shows 
that  this  contribution  is  close  in  time  to  the  specular  reflection,  a  surface  guided  wave  does 
not  have  sufficient  time  to  travel  around  the  backside  of  the  sphere  A  novel  alternative 
prompt  guided- wave  mechanism  is  from  rays  with  the  group  velocity  directed  opposite  to 
the  phase  velocity^^*^^'*  and  Fig.  12  is  generally  supportive  of  such  a  mechanism.  This 
mechanism  is  examined  below  in  Section  IV  AlO. 

7.  Ray  Synthesis  of  Leaky  Lamb  Wave  Contributions  to  Backscattering  from  Thick 
Cylindrical  Shells- 

A  ray  synthesis  as  indicated  in  this  subtitle  was  carried  out  by  Sun  and 
MarstonJ3.M5  to  test  Marston’s  earlier  analysis***®  of  the  approximate  phase  of  the 
coupling  coefficient  G/  for  leaky  Lamb  waves  on  a  cylindrical  shell.  It  is  noteworthy  that 
subsequent  analysis  by  Felsen  et  al.®  reduces  to  the  same  approximate  phase  and  amplitude 
for  G/  in  this  case  so  that  the  synthesis  can  also  be  taken  as  a  test  of  their  results.  The  ray 
picture  is  as  shown  in  Fig,  1.  Attenoon  here  is  restricted  to  the  case  of  a  stainless  steel  440 
C  shell  with  h/a  =  0,162.  Figure  14  compares  the  resulting  synthesis  (the  solid  curve)  with 
the  form  function  Ifl  from  the  exact  partial  wave  series.  The  agreement  is  good  down  to  ka 
«  7.  The  synthesis  is  quite  analogous  to  the  discussion  for  shells  given  by  Kargl  and 
Marston.**  For  future  reference  the  relevant  phase  velocity  and  dispersion  curves  are 
shown  in  Fig.  15  as  given  by  the  complex-v  plane  method. 

8.  Weakly-or-Non-Radiating  Supersonic  Lamb  Waves  on  Fluid-Loaded  Plates  and  Shells. 

Consider  again  the  conventional  leaky  ray  picture  for  backscattering  by  spheres  and 
cylinders  shown  in  Fig.  1.  The  associated  contributions  have  been  used  to  synthesize  the 
exact  f  for  thick  hollow  spherical*^*  and  cylindrical^^  shells.  In  those  ray  calculations,  the 
leaky  lamb  wave  properties  c/(kfl)  >  c  and  P/(kfl)  arc  calculated  from  the  exact  elastic 
equations  as  discussed  above  in  conjunction  with  Fig.  3.  One  curious  result  of  those 
calculations  is  that  for  certain  leaky  rays,  the  function  p/(ka)  can  exhibit  a  minimum  where 
P/  can  be  as  small  as  10*^  Np/rad  which  means  that  the  supersonic  guided  wave  only 
weakly  leaks  energy.  In  the  region  of  ka  where  P/  is  small,  the  resonances  associated  with 
that  leaky  wave  are  narrow.*  For  the  example  of  a  thick  stainless  steel  shell  with  h/a  = 

0.162  considered  in  Fig.  15,  such  a  minimum  P/  occurs  at  kn  »  18  for  the  /  =  so  wave. 
Marston*^**3.P2  showed  that  insight  into  the  cause  of  such  minima  can  be  found  by  first 


Hgare  14  -  Confbnnatioa  erf  die  ray  approximation  (solid  curve)  including  Ic^  Lamb  wave  contributions  to 
backscaticring.  The  dashed  curve  is  exact  From  Son  and  Marston  (Publication  J3). 


considering  the  particle  displacement  at  the  surface  of  a  flat  plate  in  a  vacuum.  Worltonl^ 
has  shown  that  when  the  phase  velocity  c/p  of  the  /th  symmetric  Lamb  mode  of  the  plate  is 
equal  to  the  longitudinal  wave  speed  cl  of  the  bulk  material,  then  there  is  no  normal 
displacement  at  the  plate's  surface.  For  a  supersonic  wave  to  radiate  or  leak  sound  into  the 
surrounding  inviscid  fluid,  it  is  essential  that  surface  displacements  occur  peipendicular  to 
the  plate.  Consequently,  it  is  to  be  expected  that  the  radiation  damping  vanishes  exactly 
when  c/p  -  cl  for  symmetric  waves  on  fluid  loaded  plates.  Inspection  of  the  exact 
dispersion  equations  given  in  Ref.  1 1  and  10  for  the  case  of  (inviscid)  fluids  on  one  or  both 
sides,  respectively,  shows  that  is  indeed  the  case:  when  c/p  =  cl,  there  is  no  radiation 
damping  of  the  /th  symmetric  Lamb  wave  and  c/p  does  not  depend  on  the  density  of  the 
fluid!  A  corresponding  analysis  of  the  exact  displacements  at  the  (curved)  surface  of  a  shell 
has  not  been  earned  out,  however,  it  is  to  be  anticipated  on  physical  grounds  that  locally 
normal  surface  displacements  are  essential  either  for  radiation  by  a  given  mode  or  for 
excitation  of  that  mode  by  an  incident  sound  wave.  For  leaky  waves  with  kfl  »  1  on 
shells  with  h/a  «  1,  the  effects  of  curvature  should  be  small  so  that  p/  should  be  small  for 
symmetric  waves  when  c/ «  cl-  If  curvature  is  neglected  in  Fig.  15,  the  prediction 
becomes  that  |3/(ka)  should  be  small  for  c/(ka)  such  that  c//c  «  eje  «  3.98.  Comparison  of 
Fig.  15(a)  with  15(b)  for  the  aforementioned  so  wave  at  ka  »  18  confirms  the  usefulness  of 
this  criteria.  It  is  noteworthy  that  the  existence  of  supersonic  waves  with  little  or  no 
radiation  is  contrary  to  conventional  wisdom. 

9.  Leaky  Lamb  Wave  Contributions  to  Backscattering  of  Tone  Bursts  by  Thick  Cylindrical 
Shells:  Experiments. 

During  the  early  use  of  the  6500  gallon  tank  facility,  S.  P.  Perry  (with  some 
assistance  from  D.  H.  Hughes)  carried  out  various  novel  experiments  pertaining  to  the 
backscattering  by  an  empty  thick  cylindrical  shell  for  situations  where  the  incident  tone 
burst  was  sufficiently  short  that  contributions  from  different  reflected  and  guided  waves 
were  separate  in  time.^^  The  experiment  was  analogous  to  Kargl  and  Marston's  previous 
observations  of  the  backscattering  of  tone  bursts  by  a  thick  spherical  shell^  except  that  now 
a  circular  cylinder  was  used.  The  relevant  ray  diagram  for  the  guided  wave  contributions 
studied  is  shown  in  Fig.  1.  The  empty  cylindrical  shell  used  in  the  experiments  was  made 
of  stainless  steel  304.  The  thickness-to-outer  radius  ratio  was  h/a  »  0. 17.  Figure  16(a) 
shows  the  backscattered  short  tone  burst  echo  at  a  frequency  of  160  kHz  which 
corresponds  to  ka  =  14.2.  The  backscattering  of  the  ao  leaky  Lamb  wave  is  delayed  in  time 
relative  to  the  specular  and  corresponds  to  the  earliest  or  m  =  0  circumnavigation. 


Marston's  approximation  for  IG/I  for  cylinders  ^  predicts  that  the  contribution  to  the  form 
function  for  backscattering  is 


If  I- 
ml 


J/2 

(itka) 


exp  [^-2(n  -  -  m27ip^j, 


(8) 


where  sin9/  =  c/c/  <  1.  Figure  16(b)  shows  the  predicted  Ifm/I  contributions  of  the  m  =  0 
circumnavigations  of  the  ao  and  so  waves  for  the  cylindrical  shell  under  consideration. 

Also  shown  in  the  experimental  result  obtained  from  Fig.  16(a)  by  taking  the  ratio  of  the 
guided  wave  and  specular  echo  amplitudes.  According  to  the  arrival  time,  the  guided  wave 
is  for  m  =  0  and  /  =  hq.  The  agreement  with  the  theory  supports  the  applicability  of 
Marston's  approximation  of  IG/I.  A  noteworthy  aspect  of  the  experiment  is  that  the 
frequency  was  such  that  all  Ifm/I  for  the  so  wave  are  predicted  to  be  small.  The  reason  is 
that  close  to  this  frequency  p/  is  small  for  /  =  so  as  a  consequence  of  poor  coupling 
discussed  in  Sec.  IV  A8  above.  This  results  in  a  situation  that  is  experimentally 
advantageous  since  interference  due  to  so  wave  echoes  is  predicted  to  be  weak.  In  addition 
to  this  experimental  result,  there  is  other  evidence  supporting  the  correemess  of  Eq.  (8) 
since  when  phase  infonnation  is  included  and  summations  over  the  indices  m  and  /  are 
preferred  the  resulting  guided  wave  contribution  to  the  steady  state  scattering  is  the  one 
verified  by  the  comparison  shown  in  Fig.  14?^ 

10.  High-Frequency  Scattering  Enhancement  by  Shells:  Negative-Group  Velocity  Leaky 
Wave  giving  Enhanced  Prompt  Radiation. 

Calculations  of  the  steady-state  form  function  Ifl  for  empty  spherical  shells  in  water 
having  thicknesses  h  of  less  than  approximately  10%  of  the  radius  a  display  a  prominent 
enhancement  in  a  region  at  high  frequencies.  This  enhancement  has  been  noted  in 
discussions  of  computations  by  a  few  authors  and  there  has  been  some  speculation  as  to  the 
cause.  Figure  17  gives  calculations  by  Kargl  and  Marston^>J^  for  the  case  of  an  aluminum 
shell  with  h/a  =  0.04.  The  solid  curve  shows  Ifl  which  displays  a  pronounced  peak  near  ka 
s  302.  The  dashed  curve  shows  the  specular  contribution  calculated  from  ray  theory^  for 
the  case  of  a  fluid  shell  having  the  same  thickness  and  longitudinal  sound  speed  cl  as  the 
aluminum  shell.  Both  curves  show  evidence  of  a  broad  minimum  in  the  specular 
backscattering  contribution  near  the  frequency  of  the  lowest  longitudinal  resonance  which 
is  at^Ji 
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Figure  17  --  Form  function  for  a  4%  thick  hallow  aluminum  shell  showing  a  large  high- 
frequency  enhancement  at  ka  «  302. 


BACKSCATTERING  FROM  SPHERICAL  SHELL 


Figure  18  »  Conq>uted  backscattering  of  a  20  cycie  sine  wave  bunt  with  the  carrier 
fir^uency  at  ka  ■  302,  the  large  peak  in  Fig.  17.  The  specular  reflection  at  T  ■  1  is 
followed  by  a  prompt  backscattering  enhancement  associated  wUh  the  large  peak. 


For  the  aluminum  shell  in  water,  cl  =  6.42  km/s  and  Eq.  (9)  predicts  (ka)LR  =  340.  The 
high  frequency  enhancements  predicted  in  computations  by  other  authors  are  also  at 
frequencies  close  to  but  below  the  longitudinal  resonance  condition  of  Eq.  (9). 

To  gain  insight  into  the  mechanism  for  the  enhancement,  in  Zhang's  M.S.  degree 
thesis,^^  the  backscattering  of  tone  bursts  having  carrier  frequencies  in  the  enhanced  region 
were  calculated  by  the  exact  method  of  Fourier  synthesis.  (While  the  computational  method 
is  summarized  in  Ref.  J4,  the  results  for  this  aspect  of  the  investigation  are  shown  in  Ref. 
Tl.)  Figure  18  shows  an  example  of  the  results  for  the  aluminum  shell  considered  in  Fig.  17 
when  the  carrier  ka  =  xq  is  taken  to  be  302  corresponding  to  the  peak  value  of  the 
enhancement.  The  amplitude  and  time  normalization  is  the  same  as  that  discussed  in 
conjunction  with  Fig.  5  and  there  are  20  sine  wave  cycles  in  the  incident  burst.  The  echo 
contributions  caused  by  guided  waves  that  travel  around  the  back  side  of  the  shell  do  not 
arrive  until  after  T  «  2.5.  Inspection  of  Fig.  18  (and  similar  results  for  other  burst  lengths) 
shows  that  the  high  frequency  enhancement  is  associated  with  a  build  up  beginning  by  T  » 

1.2  that  is  nearly  saturated  by  about  T  =  1.5.  Consequently  the  enhancement  is  not 
associated  with  guided  waves  that  travel  around  the  back  side  of  the  shell.  This  conclusion 
was  previously  noted  in  conjunction  with  the  discussion  of  Fig.  12. 

A  plausible  ray  diagram  that  explains  the  promptness  of  the  high  hequency 
enhancement  was  suggested  by  Marston.^  Tl^e  diagram.  Fig.  19,  should  apply  to 
situations  where  the  sign  of  the  group  velocity  of  the  guided  wave  is  negative  relative  to  the 
phase  velocity.  The  enhancement  occurs  in  the  general  frequency  region  where  a  leaky  wave 
on  the  shell  has  a  negative  group  velocity.  Such  a  reversal  of  group  velocity  was  previously 
known  for  waves  on  plates  in  a  vacuum20.2I  ^  frequency  slightly  below  that  predicted  by 
Eq.  (9).  [It  is  necessary  to  convert  to  the  units  of  Eq.  (9)  that  pertain  to  scattering 
calculations  for  a  shell  of  the  same  thickness  as  the  plate.]  The  ray  diagram  in  Fig.  19  is 
based  on  the  usual  trace  velocity  matching  condition  0  <  sin9/  =  c/c/  <  1  and  the  energy  flow 
follows  the  group  velocity  which  is  in  the  backward  direction.  The  contribution  to  the 
scattering  is  expected  to  be  prompt  since  the  energy  does  not  need  to  flow  around  the 
backside  of  the  sphere.  Indeed,  the  backside  of  the  target  is  seen  to  be  irrelevant  to  the  high 
frequency  enhancement  The  enhancement  may  be  relevant  to  the  sonar  detection  of  objects 
whose  backside  is  buried,  floating,  or  otherwise  obscured  to  sound  in  water. 

Support  for  the  negative-group  velocity  mechanism  diagramed  in  Fig.  19  was  evident 
when  a  graduate  student  working  in  this  project  D.  H.  Hughes,  discovered  a  previously 
unknown  class  of  complex-v  plane  poles  for  fluid  loaded  shells  that  become  relevant  in  the 
frequency  range  of  the  high-frequency  enhancement^*'*  Here  v  is  the  complex  partial  wave 
index  and  the  poles  correspond  to  roots  of  Dv(ka) »  0  where  Dn(ka)  =  0  is  the  denominator 
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Figure  19  -  Possible  mechanism  for  the  prompt  enhanced  backscattering  calculated  in 
Figs.  17  and  18  to  occur  for  a  narrow  range  of  ka  at  high-frequencies.  It  is  necessary  for 
there  to  exist  leaky  Lamb  waves  having  negative  group  velocities  with  only  moderate 
damping.  The  phase  velocity  of  the  wave  coupled  at  A  is  clockwise  but  the  energy  flux  is 
counterclockwise.  The  trace  velocity  condition  is  matched  for  backscattering  at  B.  From 
symmetry  there  is  superposed  ray  directed  from  B  to  A. 


Figure  20  -  ITie  mechanism  described  in  Fig.  19  may  be  facilitated  by  a  new  class  of 
conq)lex-v  plane  poles  generally  lying  in  the  4th  (or  2nd)  quadrants  of  v/ »  a/  +  ip/.  This 
figure  shows  the  loci  of  such  poles  for  a  diick  stainless  steel  shell  (widi  and  without  fluid 
loading)  as  ka  is  varied  close  to  the  si  mode  threshold  near  ka  of  (ka)LR  76.8. 


of  the  nth  term  of  the  partial  wave  series.  (See  the  discussion  of  Fig.  3  and  the  relationship 
between  roots  and  wave  properties  reviewed  in  Refs.  2, 4,  Jl,  and  J4.)  To  appreciate 
Hughes's  discovery  it  is  necessary  to  summarize  these  relationships  for  spheres.  For  the  /th 
class  of  guided  wave  at  a  frequency  ka,  the  root  of  I>v(ka)  =  0  is  denoted  by  V/  =  a/  +  iP/ 
where  the  phase  and  group  velocities  are 

cjc  =  x/(a^  +  1/2) ,  cjc  =  (da^dx)*^ ,  (10,1 1) 

where  x  =  ka.  For  the  usual  class  of  leaky  waves  considered  V/  lies  in  the  first  quadrant  of 
the  complex  plane  and  c/  and  Cg/  are  both  positive.  Hughes  discovered  a  previously 
unknown  class  of  roots  generally  lying  in  the  second  and  fourth  quadrants  of  the  complex  v 
planes.  The  root  trajectories  are  such  at  c/  and  Cg/  have  opposite  signs.  For  example,  the 
fourth  quadrant  roots  have  p/  <  0,  c/  >  0,  and  Cg/  <  0.  (Note  that  a  negative  value  of  p/  does 
not  actually  imply  a  growth  of  the  wave  with  propagation  around  the  shell  since  the  ray 
direction  is  reversed  in  Fig.  19  from  the  usual  situation  in  Fig.  1.  Since  c/  diverges  as  a/ 
approaches  •  1/2,  the  actual  mode  threshold  for  this  root  lies  outside  the  fourth  quadrant 
Poles  exhibit  a  reflection  symmetry  about  a/  =  -1/2,  p/  =  0  so  that  fourth  quadrant  poles  are 
mirrored  in  the  second  quadrant.)  Unfortunately  there  are  computational  difficulties  in 
tracking  the  new  complex  v  poles  over  the  entire  ka  interval  of  interest  for  thin  shells  though 
the  existence  of  such  poles  has  been  verified. 

For  the  purposes  of  the  present  summary,  the  numerical  results  will  be  illustrated  by 
considering  an  empty  spherical  shell  of  440C  stainless  steel  having  h/a  =  0.162  so  that  Eq. 
(9)  gives  (ka)LR  =  76.8.  This  is  the  same  shell  for  which  the  pole  trajectories  for  v  in  the 
first  quadrant  were  previously  studied  by  Ksrgl  and  Marston.2.4.5  Jl  Computations  for  this 
shell  exhibit  a  high  frequency  enhancement  anomaly  near  ka  =  72  (sec  Fig.  7  of  Ref.  Jl.). 
Figure  20  illustrates  the  relevant  root  loci  including  the  new  fourth-quadrant  root.  These 
waves  are  designated  by  /  =  si  and  /  =S2b  though  the  reader  is  cautioned  that  in  the  related 
problem  of  Lamb  waves  on  a  plate  in  a  vacuua,  the  notations  used  to  denote  the  wave  branch 
with  (c//cg/)  <  0  differ  widely  among  various  authors.  *^*22  Figure  21  shows  the  resulting 
calculation  of  c/  and  Cg/  for  the  relevant  class  of  surface  guided  waves.  Prior  to  examining 
the  trajectories  for  the  fluid  loaded  shell,  the  trajectories  are  examined  for  a  shell  in  a  vacuum 
by  letting  the  density  p  for  the  outer  fluid  (water)  vanish.  Then  the  consequences  of  taking  p 
=  water  =1.0  g/cm^  are  examined.  In  both  cases  the  frequency  units  are  taken  to  be  x  =  ka 
where  k  =  ta/c  and  c  is  fixed  at  1.479  km/s,  the  speed  of  sound  in  (real)  water.  For  the  shell 
in  a  vacuum  the  loci  of  the  si  and  S2b  trajectories  are  such  that  the  locus  of  the  c/  gives  the 
smooth  connected  curves  with  a  pitch-fork  bifurcation  at  x  =  71.6  shown  in  Fig.  21(a).  The 
corresponding  Cg/  curves,  the  solid  curves  in  Fig.  21(b),  show  the  curve  with  Cg/  >  0  with  /  = 


Si  connects  smoothly  with  curve  with  Cg/  <  0  for  /  =  S2b.  The  lower  curve  loops  back  to  Cg/ 
-  0  when  x  -  (ka)LR.  The  loop  in  tlie  group  velocity  curve  with  Cg/  <  0  is  like  the  one 
predicted  for  waves  on  a  plate  in  a  vacuum^  ^  and  the  calculations  for  shells  by  Sammelmann 
and  Hackman^  based  on  complex-ka  plane  poles  (presumably  with  negligible  fluid  loading). 
When  fluid  loading  corresponding  to  water  is  included  in  our  calculation  the  dotted  curves  in 
Fig.  21(b)  are  obtained  where  only  the  lower  curve  has  Cg/  <  0  and  (from  the  locus  in  Fig. 
20)  P/  <  0.  This  branch  corresponds  to  the  upper  continuous  curve  for  c/  in  Fig.  21(a). 
Associated  with  this  branch  (that  is  novel  for  fluid  loaded  shells)  will  be  a  backward 
propagating  wave  as  diagramed  in  Fig.  19.  Th^re  may  be  complications  in  assignment  of  the 
amplitude  for  such  ray  contributions  that  are  currently  being  investigated  though  an 
approximation  suggested  by  Marston^  appears  to  be  sufficient  for  explaining  the  order  of 
magnitude  and  general  position  of  the  enhi.  cement. 

It  is  appropriate  to  note  that  Kargl^  found  a  different  class  of  negative  group  velocity 
first-quadrant  poles  pres.'int  only  in  a  more  restricted  range  of  ka  that  for  the  present  example 
corresponds  ka  <  70.5.  There  is  computational  evid^'oee  for  thinner  sheUs  that  the  fourth 
(and  second)  quadrant  class  of  complex-v  poles  relevant  to  Fig.  20  are  the  ones  associated 
with  the  high-fi-equency  backscattering  enhaiicement. 

B .  Radiation  Mechanisms  for  Evanescent  Tone  Bursts  on  Fluid  Loaded 
Membranes  and  Plates. 

This  line  of  experiments  was  started  to  investigate  various  mechanisms  for  radiation 
by  subsonic  flexural  wave  tone  blasts  on  membranes  and  plates.  Ordinarily  the  evanescent 
wavefield  of  a  subsonic  flexural  wave  would  be  trapped  to  an  infinite  flat  plate.  In  the 
presence  of  curvature,  tunneUng  leads  to  radiation  as  diagramed  in  Fig.  2  for  a  guided  wave 
that  is  only  slightly  subsonic  near  the  coincidence  frequency.  Radiation  is  possible  for  a 
flat  plate  or  membrane  in  response  to  a  perturbation  of  the  wavefield.  This  is  diagramed  in 
Fig.  22  for  the  case  of  an  external  obstruction.^^  Radiation  into  the  fluid  may  similarly 
result  when  a  subsonic  wave  impinges  on  internal  structures  such  as  a  rib  or  at  a  sudden 
transition  in  fluid  loading  such  as  when  a  wave  on  a  plate  crosses  an  interface  into  water. 

Since  the  problems  of  interest  concerns  subsonic  flexural  waves,  it  was  decided  to 
excite  such  waves  directly  with  an  electromagnetic  acoustic  wave  transducer  (EMAT)  rather 
lhan  by  driving  the  reciprocal  of  the  process  to  be  studied  with  an  external  sound  wave. 

When  such  transducers  are  used  for  surface  waves,  the  typical  application  is  for 
frequencies  above  ~  1  MHz.23  In  the  present  application  the  typical  frequency  is  ~  100 
kHz  or  less  and  the  width  of  the  region  to  be  excited  exceeds  ~  5  inches.  The  initial 


experiments  for  this  aspect  of  Marston’s  research  program  involved  the  proof  of  principal 
for  the  use  of  meanderline  EMATs  to  excite  low  frequency  flexural  waves.  The 
experiments  were  done  on  a  stretched  aluminized  Mylar  membrane  in  air.^®  jhe 
resulting  tone  bursts  driven  on  the  membrane  were  then  used  to  develop  wave-number 
selective  optical  and  electrical  displacement  transducers  for  detecting  such  bursts.^^  Some 
of  those  transducers  were  subsequently  found  useful  for  detecting  flexural  wave  tone 
bursts  on  plates. 

Flexural  tone  bursts  on  aluminized  Mylai'  that  were  subsonic  relative  to  the 
surrounding  air  were  used  to  conHim  one  of  the  radiation  processes  illustrated  in  Fig.  22. 

The  tone  bursts  had  a  frequency  f  =  4  kHz,  a  wavelength  along  the  Mylar  X/  =  1.5  cm,  and 
a  phase  velocity  c/  =  62  m/s.  The  predicted  evanescent  pressure  field  near  the  membrane 
is2432J9 

p  =  Pq  exp(-kyy)  exp(ik^x  -  icDt) ,  ky  =  (k^  -  (oVc))^^  (12,13) 

where  here  c  =  340  m/s  becomes  the  speed  of  sound  in  air,  kx  =  27cA/,  and  ky  =  l/(2.43 
mm).  The  problem  of  diffraction  of  an  evanescent  sound  wave  by  a  (rigid)  sharp  edge  has 
been  considered  for  the  corresponding  optics  problem  by  Felsen^^  and  others.  The 
analysis  should  apply  to  the  present  situation  when  the  wave  scattered  by  the  edge  is  only 
weakly  rescatiered  by  the  membrane.  The  two-dimensional  result  may  be  put  in  the  form 

'Pscat'  ”  Pq  .  (14) 

where  h  is  the  distance  of  the  edge  from  the  membrane,  r  is  the  distance  of  the  field  point  from 
the  edge,  and  A  is  a  constant.  To  test  this  result,  a  shaip  metal  blade  was  used  as  a  barrier  and 
a  probe  microphone  was  mounted  a  small  distance  away  as  shown  in  Fig.  23(b).  Tone  bursts 
were  used  that  typically  exceeded  10  cycles  in  length  but  were  sufficiently  short  to 
discriminate  in  time  certain  background  noise  sources.  (In  addition,  certain  repeatable 
background  noises  were  subtracted  from  all  of  the  temporal  records.)  The  blade  distance  from 
the  membrane  h  was  varied  while  the  microphone  was  held  fixed.  It  was  not  possible  to 
achieve  acceptable  signal-to-noise  levels  with  the  microphone  in  the  true  farfield  of  the  barrier 
(car/c  »  1).  The  mi*  ophone  was  sufficiently  distant,  however,  that  variations  in  r  had  a 
negligible  effect  on  the  received  signal  for  a  gap  h  up  to  about  1  mm.  The  observations.  Fig. 
23(a),  confirm  that  the  scattered  signal  is  largest  when  the  blade  is  close  to  the  membrane  but 
not  in  contact  with  it.  Furthermore,  the  signal  decreases  with  increasing  h  roughly  in 
proportion  to  exp(-kyh)  as  predicted  from  a  theory  that  neglects  rescattering  by  the  thin 
membrane.  The  experiment  coniums  that  perturbations  of  an  evanescent  acoustic  wavefield 
do  indeed  result  in  scattering.I9J2,M15 


41 


Figure  22  '•  Radiation  to  the  farfield 
can  occw  if  objects  are  placed  in  die 
evanescent  wavefield  of  a  subsonic 
wave  packet  (from  Marston, 
Publication  62). 
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Figure  23  --  Confirmation  of  the 
dimaction  process  illustrated  in  Fig. 
22  but  for  a  subsonic  wave  packet  on 
a  membrane  in  air. 
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The  principal  recent  effort  has  been  to  establish  a  technical  capability  of  exciting  tone 
bursts  of  waves  on  metal  plates  near  and  below  the  coincidence  frequency  with  water  on  both 
sides  of  the  plate  and  to  probe  the  pressure  field  near  the  plate.  Some  of  the  motivations  for 
establishing  this  capability  will  be  evident  from  the  discussion  below  and  from  the  research 
proposed  in  Sec.  m.  The  approach  has  three  technical  aspects:  (a)  water  tank  facility,  (b) 
nods  selective  EMAT  source,  and  (c)  probes  of  the  acoustics  field  in  the  water.  Progress  is 
summarized  as  follows:  (a)  Tank:  to  avoid  interference  with  scattering  experiments  carried 
out  in  the  large  12  ft  diameter  x  8  ft  deep  tank  mentioned  earlier,  a  second  wooden  tank,  8  ft 
diameter  by  8  ft  deep  was  assembled  adjacent  to  the  large  one.  The  tanks  use  a  common  metal 
platform  to  facilitate  access  to  the  experimental  region  and  each  is  mechanically  isolated  from 
adjacent  structures.  See  Section  V.  (b)  EMAT  Sources:  The  source  was  developed  in 
cooperation  with  Dr.  Bruce  W.  Maxfield  (Industrial  Sensors  and  Actuators,  San  Leandro, 

CA)  drawing  upon  his  previous  experience.23  A  functional  diagram  is  shown  in  Fig.  24. 

Cfhe  transducer  is  mounted  above  the  water  and  the  plate  bends  down  into  the  water  in  a 
manner  to  be  described  subsequently.)  A  tone  burst  of  current  is  applied  to  antisymmetric 
coils  adjacent  to  the  plate  but  not  in  contact  with  it.  (The  actual  coils  used  differ  in  detail  from 
the  meander  coils  diagramed.)  The  plate  thickness  (1/8  inch)  is  much  thicker  than  the  skin 
depth  for  the  typical  frequencies  used  (~  30  kHz  and  ~  100  kHz).  The  coils  induce  opposing 
currents  on  each  side  of  the  plate  such  that  current  loops  are  completed  near  the  edges  of  the 
plate.  Those  cunents  result  in  oppositely  directed  Lorentz  forces,  which  if  pr  operly  balanced 
produce  a  pure  couple  (or  torque)  that  cause  the  plate  to  bend.  The  coils  are  periodic  with  a 
period  2b  =  1.5  cm.  Consequently,  counter  propagating  antisymmetric  Lamb  waves  are 
excited  having  a  wavelength  of  2b.  (Wavelengths  that  are  multiples  or  certain  rational 
fractions  of  2b  have  also  been  excited.)  When  the  frequency  is  selected  so  that  2b  matches  the 
wavelength  of  a  mode,  the  waves  are  reinforced  by  adjacent  coils  such  that  the  amplitude  of 
the  radiated  wave  is  increased  by  the  number  N  =  8  of  complete  coil  cycles.  The  magnetic 
field  strength  B  is  4  kGauss  and  with  the  available  current  source,  the  displacement  amplitude 
of  the  radiated  Lamb  wave  tone  burst  is  about  10"^  m.  The  amplitude  was  inferred  by  using  a 
differential  wavenumber  selective  capacitor  microphone  as  originally  proposed,  (c)  Probes  of 
the  Acoustic  Field:  The  wave  packets  propagate  from  the  plate  in  air  to  the  portion  of  the  plate 
in  water.  This  can  result  in  a  superposition  of  leaky  and  evanescent  waveflelds  near  the  plate. 
Hydrophones  are  used  to  probe  this  wavefield  which  are  positioned  by  a  two-dimensional 
stepper  motor  positioning  system.  Figure  36(b)  of  Section  V  is  a  photograph  of  the  EMAT 
and  plate  and  positioning  system  as  mounted  above  the  tank. 

Figure  25  illustrates  one  of  the  experimental  coiiEgurations  studied  for  a  flexural 
wave  on  a  0.125  inch  thick  aluminum  plate.  (By  setting  the  phase  velocity  in  a  vacuum  of 


the  ao  Lamb  wave  to  the  speed  of  sound  c  in  water,  the  coincidence  frequency  fc  is 
estimated  to  be  89  kHz.)  Figure  26(a),  (b),  and  (c)  illustrates  representative  time  records 
with  the  hydrophone  located  respectively  in  the  positions  labeled  1, 2,  and  3  in  Fig.  25. 
The  hydrophone  separation  from  the  plate  for  positions  1  and  2  was  about  3  mm.  The 
records  are  for  12  cycles  of  drive  to  the  EMAT  with  a  carrier  frequency  of  32  kHz  where 
the  phase  velocity  c/  is  subsonic  with  respect  to  the  velocity  of  sound  c  in  water.  The  early 
time  signal  on  the  left  is  due  to  direct  electrical  pickup  of  current  tone  burst  that  drives  the 
EMAT.  The  larger  of  the  bursts  near  the  center  and  right  of  (a)  and  (b),  respectively,  are 
primarily  from  pressure  variations  in  water  resulting  from  the  passage  of  a  flexural  wave 
packet  on  the  fluid-loaded  plate.  Its  group  velocity  Cg  is  1 .6  mnV)is  as  determined  from  the 
time  shift  and  difference  in  propagation  distances  between  (a)  and  (b)  and  agrees  with 
prediedons.  The  signal  near  the  left  side  of  (c)  has  a  delay  time  characteristic  of 
propagation  directly  from  where  the  plate  enters  the  water,  E  in  Fig.  25,  to  the  hydrophone 
at  position  3.  (The  calculation  of  the  arrival  time  includes  the  delay  for  the  wave  packet  to 
propagate  from  the  EMAT  to  point  E.)  The  decay  of  amplitude  with  the  distance  y  from  the 
plate  is  weaker  than  predicted  by  Eq.  (12)  which  is  to  be  expected  if  the  signal  is  radiated 
from  a  line  source  at  E  instead  of  decaying  as  an  evanescent  wave.  This  radiation 
mechanism  is  described  below.  It  is  remarkable  that  the  signal  in  Fig.  26,  as  well  as  those 
for  Fig.  27  described  below,  confirm  that  the  properly  aligned  EMAT  is  highly  selective  in 
the  plate  mode  excited  so  that  only  the  ao  mode  is  driven. 

Figure  27  displays  records  for  excitation  of  the  EMAT  with  16  cycles  at  106  kHz. 
The  hydrophone  positions  for  (a)  and  (b)  are  positions  1  and  2  of  Fig.  25  and  the  wave 
packets  again  are  primarily  from  pressure  variations  in  water  resulting  from  the  passage  of 
a  flexural  wave  packet  on  the  fluid-loaded  plate.  The  wave  packet  in  (c)  corresponds  to 
direct  radiation  form  a  line  source  at  E  but  the  hydrophone  position  is  not  shown  in  Fig.  25 
being  at  a  depth  of  50  cm  and  a  horizontal  distance  from  the  plate  y  of  50  cm. 

From  the  delay  times  for  the  wave  packets  shown  in  Figs.  26  and  27  (and  other 
records  not  shown)  the  following  is  evident:  when  the  hydrophone  is  far  from  the  plate, 
the  signal  is  dominated  by  direct  radiation  from  where  the  plate  enters  the  water.  We  call 
this  radiation  transition  radiation  because  of  a  similar  terminology  used  to  describe  the 
electromagnetic  radiation  from  relativistic  charged  particles  when  they  cross  the  surface  of  a 
dielectric.  A  detailed  inspection  of  the  records  in  (a)  and  (b)  of  Figs.  26  and  27  suggests 
that  those  records  also  contain  a  transition  radiation  component  That  the  wavefield 
consists  of  a  superposition  is  in  agreement  with  the  following  considerations. 

First  consider  the  wavefield  near  an  infinite  plate  for  which  the  pressure  vanishes 
on  an  unbounded  peipendicular  plane  corresponding  to  the  free  water  surface  in  Fig.  25. 


Amplflude  (v)  Amplrtude  (v)  Amplllude  (v) 


The  wavefield  that  meets  those  boundaiy  conditions  in  the  presence  of  a  plate  corresponds 
to  that  of  a  harmonic  moment  (or  torque)  applied  at  the  line  where  the  plate  crosses  the  free 
surface  of  the  water.  The  wavefield  given  by  differentiation  (with  respect  to  the  source 
coordinate)  of  the  Green's  function  for  a  line  source  on  a  fluid-loaded  plate  as  described 
e.g.  by  Feit^  but  modified  to  include  fluid  on  both  sides  of  the  plate.  It  is  evident  from 
Feit's  analysis  for  excitation  below  the  true  coincidence  firquency  that  the  wavefield  will 
consist  of  a  superposition  of  an  evanescent  wavefield  with  that  of  a  wave  propagating  in 
water  from  the  line  source.  The  direct  wave  results  from  nonpropagating  bending  of  the 
plate  localized  near  the  line  source.  In  the  present  case,  localized  bending  at  the  free 
surfaces  give  rise  to  the  transition  radiation.  A  quantitative  understanding  of  wavefield 
superposition  in  the  water  is  needed  before  meaningful  scattering  experiments  of  the  type 
illustrated  in  Fig.  22  can  be  completed.  This  wavefield  appears  to  be  a  canonical  problem 
of  general  interest:  the  radiation  from  a  line  moment.  The  experiment  has  the  added 
complication  of  the  finite  width  of  the  plate. 

Summary  of  experimental  results:  The  diffraction  of  evanescent  waves  by  a 
barrier  in  air  was  demonstrated  and  the  variation  of  signal  level  with  barrier  position  is  in 
agreement  with  a  simple  approximation  that  neglects  the  rescattering  of  the  waves  by  the 
membrane.  The  experiments  for  flexural  waves  on  a  plate  in  water  are  more  complicated 
due  to  the  nature  of  the  wavefield  generated  when  the  flexural  wave  packet  propagates 
along  the  plate  into  the  water.  The  EMAT  driver  for  the  plate  waves  was  found  to  be 
highly-mode  selective  and  this  greatly  simplifies  the  interpretation  of  the  experiments.  In 
this  case  the  investigation  of  the  wavefield  in  the  water  even  without  any  scatterer  present  is 
of  importance  due  to  the  contribution  by  nonpropagating  bending  of  the  plate  at  the  free 
surface  of  the  water. 


C.  Acoustical  Four-Wave  Mixing  in  Suspensions 

For  his  Ph.D.  dissertation  completed  in  1992,  H.  J.  Simpson  investigated  a  novel 
mechanism  for  the  interaction  of  sound  with  sound.^  The  new  coupling  mechanism  relies 
on  the  slow  response  to  the  acoustic  radiation  pressure  of  pump  waves  by  the  acoustic 
medium,  a  suspension  of  particles  in  water.  Three  acoustic  waves,  two  pump  wave  and  a 
higher  frequency  probe  wave,  interact  to  produce  a  fourth  wave.  The  coupling  mechanism 
is  analogous  to  optical  four-wave  mixing  in  photo-refractive  materials.^^*^?  7^5  summary 
of  results  given  below  illustrates  some  of  the  nonlinear  properties  of  this  system,  however 
the  reader  is  referred  to  the  dissertation  for  technical  and  theoretical  details.  The 


suspensions  studied  were  dilute,  the  volume  fraction  of  suspended  particles  being  typically 
1.6  X  lO'^.  Nevertheless,  the  wave  generated  by  the  interaction  was  not  weak  and  it  could 
be  easily  observed  with  a  large  signal-to-noise  ratio  for  a  properly  aligned  experiment. 

Figure  28  is  a  diagram  of  the  apparatus.  The  pump  wave  transducers  are  driven 
continuously  and  have  a  fixed  frequency  of  typically  800  kHz  corresponding  to  a  resonance 
for  multiply  reflected  waves  within  the  chamber.  [The  losses  are  sufficiently  small  (since 
the  transducers  are  backed  by  air)  that  it  is  only  necessary  to  drive  the  upper  transducer  in 
the  flgure  such  that  the  second  pump  wave  is  produced  by  reflections  from  the  lower 
trans^lucer.  The  lower  pump  transducer  is  used  to  monitor  the  pump  level.]  The 
suspension  consist  of  25  pm  diameter  polymer  (Latex)  spheres  in  a  saline  solution  such 
that  the  spheres  were  neutrally  buoyant  In  response  to  the  radiation  pressure  of  the 
standing  wave  some  of  the  spheres  migrate  toward  pressure  modes  of  the  standing  wave. 
The  amount  of  enrichment  of  the  local  number  density  of  spheres  is  limited  by  their 
Brownian  motion.  A  third  wave  is  introduced  in  the  form  of  a  high  frequency  (2  -  10 
MHz)  probe  wave.  An  outgoing  fourth  wave  is  produced  at  the  same  frequency  as  tlie 
probe  wave  when  the  following  Bragg  condition  is  met 

kAcos0  =  2jcm,  m=  1,2,3...  (15) 

where  0  is  the  angle  of  incidence  (for  the  probe  wave)  and  reflection  (for  the  outgoing 
wave),  A  is  the  wavelength  of  the  of  the  pump  wave  such  that  d  =  A/2  =  0.98  mm  is  the 
separation  between  enriched  layers,  and  k  =  2jtf/co  =  27t/A.probe  is  the  wavenumber  of  the 
probe  wave.  The  suspension  is  confined  by  acoustically  transparent  walls  of  Mylar  such 
that  the  pump  and  probe  transducers  both  of  water .  The  apparatus  is  aligned  such  that  0  = 
4r.  The  probe  signal  consists  of  gated  tone  bursts  that  are  sufficiently  long  to  simulate 
steady-state  Bragg  scattering  but  they  are  short  enough  for  the  radiation  pressure  of  the 
probe  wave  to  have  no  significant  effect.  The  absolute  amplitude  reflection  coefficient  IRI 
was  obtained  with  an  appropriate  calibration  achieved  by  rotating  the  receiver  transducer  so 
as  to  directly  view  the  transmitted  probe  wave.  Figure  29  gives  a  diagram  of  the  principal 
systems  used  in  the  experiment. 

Figure  30  shows  the  received  tone  bursts  corresponding  to  Bragg  orders  at  probe 
frequencies  f  of  3.27  and  4.35  MHz.  For  comparison  a  weaker  3.77  MHz  received  signal 
is  shown  that  is  representative  of  the  situation  where  the  Bragg  condition  is  not  satisfied. 
Figure  31  shows  a  response  surface  obtained  by  scanning  f  and  measuring  IRI.  The  peak 
amplitude  P  of  the  pump  standing  wave  is  also  varied  beginning  with  P  =  162  Pa  and 
increasing  to  9350  Pa.  The  Bragg  peaks  are  manifested  clearly  as  the  parallel  rows  on  this 
surface. 


BRAGG  SCATTERING  FROM  FOUR  WAVE  MIXING 
IN  NONLINEAR  MEDIA  OF  25  \m  LATEX  PARTICLES 
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Figure  28  --  A  schematic  of  the  experimental  setup  used  to  study  Bragg  scattering 
from  a  grating  of  latex  particles.  The  cell  contains  a  7%  Saline  solution  used  to  make  the 
25  |un  latex  particles  neutrally  buoyant,  and  is  enclosed  on  the  sides  of  the  cell  with  2  mil 
Mylar  to  allow  probing  of  the  fonned  grating  with  die  PVDF  transducer.  The  scattered 
wavefield  is  then  scanned  widi  the  Lithium  Sulfate  hydrophone. 
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The  development  of  a  theory  which  was  capable  of  explaining  quantitatively  the 
observed  values  of  IRI  will  now  be  summarized.^  Initially  the  case  of  an  equilibrium 
grating  will  be  examined  for  a  fixed  pump  amplitude  P  and  the  situation  of  an  unbounded 
grating.  As  a  consequence  of  the  small  volume  fraction  of  spheres  the  grating  could  be 
treated  as  collection  of  noninteracting  particles  in  a  spatially  oscillating  potential  well  U(z) 
resulting  from  the  single  particle  radiation  pressure  of  the  standing  wave.  The  local  particle 
number  density  np(z)  =  A  exp(-  U(z)/kBT)  where  ks  is  Boltzman's  constant,  T  is  the 
absolute  temperature,  and  A  is  a  constant  coefficient  determined  from  conservation  of  the 
total  number  of  spheres.  The  coordinate  z  is  parallel  to  the  pump  wave  vector  and 
perpendicular  to  each  pressure  nodal  surface  of  the  standing  wave.  Then  the  effective 
variation  of  sound  speed  c(z)  is  calculated  using  a  standard  mixture  rule  that  simplifies  in 
the  present  case  where  the  mass  density  of  the  particles  and  host  liquids  are  the  same.  The 
problem  becomes  one  of  determining  the  Bragg  scattering  from  a  specified  periodic 
variation  c(z).  Two  approaches  are  used. 

In  the  frrst  approach  the  local  variation  in  c(z)  is  approximated  as  the  profile  of  an 
Epstein  layer  for  each  spatial  cycle.  That  approach  has  the  advantage  that  the  reflection 
coefficient  Rr  for  such  a  single  layer  is  known  at  all  frequencies.  It  has  the  limitation  that 
the  Epstein  layer  best  approximates  the  calculated  local  c(z)  oscillations  at  large  pump 
amplitudes.  An  additional  approximation  is  made  that  is  analogous  to  a  Bom 
approximation  in  quantum  scattering  theory.  This  is  to  take  the  probe  beam  amplitude  as 
uniform  throughout  the  medium.  (Scattering  of  reflected  waves  off  of  other  layers  is  also 


neglected.)  The  resulting  multiple  layer  reflection  coefficient  is  approximated  as 

sin^N  kdcos9^ 

,R(.,9),  =  VN,  R^(k.9, 


where  N/  is  the  total  number  of  reflecting  layers  crossed  by  the  incident  probe  beam  and  v 


«•  1/3  is  an  q)erture  coefficient  that  is  calculated  geometrically  from  the  volume  of  overlap 
between  the  probe  beam  and  the  cylindrical  volume  viewed  by  the  probe  receiver.^ 

In  the  second  approach  the  acoustical  refractive  index  [c(yc(z)],  where  cq  is 
constant,  is  expressed  by  a  spatial  Fourier  series.  Then  analytical  results  from  optics^^  are 
used  for  Bragg  scattering  by  a  volume  grating  with  specified  Fourier  coefficients.  It  is 
necessary  for  the  Bragg  condition,  Eq.  (15),  to  be  satisfied. 

Prior  to  conqiarison  with  experimental  results  it  is  appropriate  to  con^are  the  two 
theoretical  approaches  by  restricting  attention  to  frequencies  f  that  meet  the  Bragg 
condition,  Eq.  (15).  Figure  32  shows  that  the  Epstein  layer  approximation  (dashed  curve) 
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The  General  Comnonent  Conflguration  for  Acoustical  FWM 
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Figure  29  -- Component  used  to  generate,  control  and  measure  bulk  four- wave 
mixing.  An  XT  Qone  was  used  to  control  the  rest  of  the  equipment  via  IEEE488  and 
RS232  interfaces.  Two  sets  of  function  generators  and  power  amplifiers  were  used  to 
drive  the  pump  and  probe  transducers,  and  various  low  noise  signal  processing  equipment 
was  used  to  measure  the  standing  wavefield  and  scattered  probe  wavefield. 


Received  Scattered  Signal  Off  Formed  Grating 


Received  Tone  Burst  at  3.27  MHz  Off  Grating 


Figure  30  —  Sampled  scattered  tone  burst  at  various  frequencies  on  and  off  the 
Bragg  condition.  The  RMS  voltage  of  the  viewed  dme  trace  is  obtain  using  an  algorithm 
within  the  digital  oscilloscope.  The  same  algorithm  is  used  to  measure  the  pump  voltages. 
The  top  time  trace  shows  the  scattered  response  for  a  firequency  that  satisHes  the  Bragg 
condition.  Thus  for  the  geometry  of  the  system  and  the  frequency  of  the  probe  wavefield, 
constructive  interference  occurs  for  the  reflection  firom  the  ensonified  layers.  For  a  slightly 
higher  fiequency,  as  shown  in  the  middle  time  trace,  the  Bragg  condition  is  no  longer 
satisfied.  The  bottom  time  trace  shows  the  next  frequency  that  satisfies  the  Bragg 
condition. 


53 

Experimental  Normalized  ReHection  Coefficient 
as  Frequency  and  Pressure  are  Varied 
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Figure  31  —  Experimental  reflection  coefficient  plotted  as  the  pun^)  amplitude 
and  the  probe  fiequcncy  are  varied.  The  punp  pressure  is  incremented  from  162  Pa  to 
9330  Pa  by  102  Pa,  where  the  pressure  is  maintained  constant  for  2  minutes.  After  the  2 
minutes,  the  frequency  response  of  the  grating  is  measured  and  the  pump  pressure  is  then 
incremented  tot  he  next  value. 


Figure  32  --  Theoretical  reflectivity  IRI  at  the  4th  Bragg  order  for  an  equilibrium  particle 
suspension  plotted  as  a  function  of  pump  pressure  P.  The  dashed  curve  approximates  each 
pressure  node  as  an  Epstein  layer  while  the  solid  curve  use  a  more  exact  Fourier  method. 
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Figure  33  ~  The  solid  curve  is  the  measured  reflectivity  of  die  susi^nsion  as  a  function  of 
the  probe  firequency  f.  The  dashed  curve  is  the  theoretical  reflectivity  using  the  Epstein 
layer  approximation  for  the  pun^  pressure  P  «  12200  Pa. 


is  remarkably  close  to  the  presumably  more  exact  method  grating  based  on  the  Fourier 
coefFicients.29  The  plot  corresponds  to  the  m  =  4  Bragg  peak  at  f  =  4.44  MHz  for 
conditions  representative  of  the  experiment  At  low  pump  amplitudes  P  the  grating 
reflectance  IRI  is  small  because  the  variations  in  IU(z)l  are  «  ksT  and  the  migration  of 
particles  to  pressure  nodes  is  weak.  As  P  is  increased,  IRI  grows  and  saturates  at  close  to 
the  same  value  for  both  theories.  Saturation  is  associated  with  the  condition  lAUI »  kfiT 
such  that  the  thermal  spread  of  particles  along  each  pressure  node  is  small.’^ 
Unfortunately,  it  is  not  possible  to  observe  in  equilibrium  the  rise  of  IRI  below 
approximately  5000  Pa  because  the  time  for  particles  to  migrate  to  the  pressure  nodes  is 
prohibitively  long  due  to  the  viscous  drag  on  the  particles.  Thus  neither  curve  of  Fig.  32 
actually  reproduces  the  corresponding  fixed  f  slice  through  Fig.  31  at  low  amplitudes 
where  the  grating  has  not  reached  equilibrium. 

For  the  reasons  noted  above  a  quantitative  test  of  the  equilibrium  theory  for  the 
conditions  of  our  apparatus  must  be  carried  out  at  large  pump  amplitudes  where  IRI  has 
nearly  saturated  and  the  different  approximations  are  in  good  agreement  Consequently,  it 
is  sufficient  to  examine  the  dependence  of  IRI  on  the  probe  frequency  f  at  fixed  angle  6  and 
pump  pressure  P.  Figure  33  shows  a  representative  example  of  such  a  comparison 
between  experiment  (solid  curve)  and  the  approximate  theory  (dashed  curve)  given  by  IRI 
from  Eq.  (16).  No  adjustable  scaling  parameters  were  used  in  the  comparison.  While  the 
measured  peak  amplitudes  are  slightly  below  the  predicted  amplitudes,  the  width  of  the 
observed  Bragg  lines  are  slightly  wider  than  the  predicted  ones.  This  may  be  a 
consequence  of  a  plausible  lack  of  perfect  spatial  periodicity  assumed  in  the  theory.  It  is 
evident  by  inspection  that  the  areas  of  the  measured  and  calculated  Bragg  lines  are 
similar. 'Fhe  approximate  linearity  of  the  predicted  peak  amplitude  with  probe  riequency 
f  is  a  result  of  the  single  Epstein  layer  reflectivity  IRrI  in  Eq.  (16)  being  linear  in  f  for  a  thin 
layer.T2  When  the  Bragg  condition  is  met,  Eq.  (16)  shows  that  IRI  exceeds  IRrI  by  a  factor 
vN/ »  1  and  from  Fig.  33  is  seen  to  be  as  large  as  0.01. 

A  theoi>’  for  the  time  dependence  of  the  grating  reflectivity  during  initial  grating 
formation  was  also  developed  and  confirmed  witli  experin>ents.T2  The  situation  considered 
is  the  initial  build  up  of  the  grating  from  an  homogeneous  suspension  when  the  pump 
amplitude  P  is  stepped  from  zero  to  a  constant  value.  The  theory  accounts  for  the  initial 
viscosity-limited  migration  of  spheres  to  the  antinode  but  not  the  q>proach  to  equilibrium. 
The  scattering  theory  is  based  on  a  time-dependent  expansion  of  the  spatial  Fourier 
coefficients  and  (like  the  solid  curve  in  Fig.  32)  is  restricted  to  probe  frequencies  f  meeting 
the  Bragg  condition.  Figure  34  show  the  comparison  between  the  theory  (smooth  curves) 
and  the  experiment  (noisy  curves)  for  different  values  of  P.  The  larger  P,  the  more  rapidly 


Pump  Pressure  Variation  ivith  Time 


Measured  and  Calculated  Reflectivity  Response  to  Pump  Field 


Figure  34  -- Sudden  turn  on  ofpump  pressure  field  and  measured  diffraction 
efficiency  of  the  forming  grating.  The  curves  are  the  measured  and  calculated  reflection 
coefficients  for  the  transient  grating  at  times  after  the  introduction  of  a  standing 
wavefield.  Experimentally,  the  scattered  amplitude  is  measured  every  3  seconds  and  the 
reflection  coefficient  is  calculated.  The  pump  pressure  amplitudes  used  to  calculate  the 
reflection  coefficients  arc  21500  Pa,  19400  Pa,  16900  Pa,  15000  Pa,  12800  Pa  and  10700 
Pa  in  order  from  left  to  right. 


IRI  increases  with  time.  The  experiments  were  limited  to  the  fourth  Bragg  order  (m  =  4  in 
Eq.  (15))  that  for  the  condition  of  this  experiment  gave  a  probe  frequency  f  =  4.24  MHz. 
Inspection  of  Fig.  34  suggests  a  general  but  not  perfect  agreement  between  the  experiment 
and  the  theory.  The  experiment  demonstrates  the  ability  to  measure  much  weaker  four- 
wave  mixing  processes  than  for  the  example  shown  in  Fig.  33. 

In  summary,  a  new  class  of  acoustical  four-wave  mixing  process  has  been 
demonstrated  along  with  approximate  theories  for  the  signal  strength.  The  process 
considered  is  nonlinear  in  tliat  if  either  pump  wave  is  removed  there  is  no  signal  as 
indicated  by  the  negligible  observed  signal  at  zero  time  in  Fig.  34.  A  potential  application 
of  this  process  is  an  all  acoustical  probe  of  suspensions  since  the  steady-state  (Fig.  33)  and 
time  dependent  responses  (Fig.  34)  depend  on  different  suspension  properties.  In  addition 
to  particle  size,  host  viscosity,  and  volume  fraction,  the  acoustics  signal  may  prove  useful 
for  probing  specific  condensed  matter  processes.  For  example,  if  the  volume  fraction  was 
increased  by  a  factor  of  150  from  the  conditions  observed,  the  interaction  between  particles 
should  affect  the  sicnal  since  it  would  no  longer  be  possible  for  the  particles  to  migrate  to 
layer  along  the  pressure  nodes  without  coming  into  contact.^  For  suspended  particles 
much  smaller  tlian  tiie  wavelengths  of  frie  pump  and  probe  waves,  the  Bragg  scattering 
should  depend  onh  w-eakly  on  theTatracie  shape.  Other  media  that  should  exhibit 
acoustical  four-wave  mixing  inci.uaE  emulsi^;;  and  microemulsions.^l 

It  is  notewortliy  that  after  these  experiments  were  started  an  experiment  was 
published  by  Sato  et  al.32  that  also  made  use  of  the  response  of  a  suspension  to  radiation 
lJich!:iiit  III  (/kkIucc  a  scattered  acoustic  wave.  That  experiment  had  a  more  complicated 
geometry  than  ours  (Fig.  28  and  29)  and  it  made  use  of  additional  "fixing"  standing 
acoustic  wave  to  inhibit  particle  migration  during  readout  by  the  probe  wave.  Another 
difference  i  s  tint  the  pump  waves  were  off  during  readout.  Simpson's  experiments  and 
analysis  clarifies  the  physical  mechanisms  and  interpretation^  of  the  experiments  by  Sato 
et  al. 

Chapter  6  ('f  Simpson's  dissertation  describes  experimental  and  analytical  results 
pertaining  to  a  novel  mechanism  for  the  three-wave  interaction  of  sound-with-sound  at  a 
free  surface.  That  work  was  carried  out  early  in  Simpson's  research  work  with  the  support 
of  ONR  Grant  N0(X)14-85-C-0141.  It  follows  from  an  earlier  line  of  research  on  the 
properties  of  a  localized  nonlinear  layer. 


D.  Supplemental  Scattering  Research 


In  addition  to  the  principal  topics  of  scattering  research  described  in  Sections  A  and 
B,  supplemental  topics  were  examined  with  the  partial  support  of  this  grant.  The 
investigations  mostly  concerned  examples  of  intermediate  difficulty  that  served  to  reinforce 
approaches  taken  on  primary  research  objectives.  Most  cases  involve  extensions  of 
research  previously  supported  by  ONR  Physics  Acoustics  Core  Program  Contract 
N00014-8S-C-0141.  Perhaps  the  most  important  accomplishment  in  this  category  was  the 
completion  of  a  survey  of  geometrically  based  methods  in  scattering  emphasizing  the  shape 
of  wavefronts.®^  Figure  35  shows  the  table  of  contents. 

1 .  Reflection  by  Fluid  Shells  and  Novel  Background  for  Elastic  Shells. 

Previous  research  by  Kargl  and  Marston,^-^  gave  a  novel  ray  representation  for  the 
contribution  to  the  backscatteiing  due  to  the  multiple  internal  reflections  within  an  empty 
spherical  fluid  shell.  The  result  was  found  to  be  in  good  agreement  with  novel  exact 
calculations  based  on  a  partial-wave  series  (PWS).  At  low  ka  the  phase  of  the  scattering  in 
similar  to  that  of  a  perfectly  soft  sphere  while  for  large  ka  (except  near  longitudinal 
resonances,  Eq.  (9)),  the  phase  becomes  that  of  a  rigid  sphere.  Consequently  the 
theoretical  case  of  a  vacuum  filled  spherical  shell  is  a  candidate  for  a  general  "background" 
contribution  to  "isolate"  the  elastic  response  of  shells.  An  investigation  was  carried  out  by 
Zhang  of  how  this  background  serves  to  isolate  the  prompt  contribution  to  the  high- 
frequency  enhancement  of  an  elastic  shell  for  the  case  of  tone  burst.”^*  Computations 
(Appendix  E  of  Ref.  T1 .)  demonstrate  that  with  this  background  the  prompt  elastic 
contribution  appears  to  be  isolated  from  other  contributions  for  the  specific  case  of  a  4% 
thick  aluminum  shell  considered.  Zhang^l  also  found  from  computations  of  the 
backscatteiing  of  tone  bursts  by  fluid  shells  that  certain  features  for  ka  close  to  the 
longitudinal  resonance,  Eq.  (9),  were  well  explained  by  the  ray  model  of  Kargl  and 
Marston.5  In  related  work,  Marston®  extended  the  previous  ray  analysis^  to 

multiple  internal  reflections  within  an  empty  fluid  cylinder.  With  assistance  from  Sun, 
comparisons  were  made  with  the  exact  backscatteiing  as  given  by  a  PWS  for  the  empty 
fluid  cylinder.  The  level  of  agreement  is  similar  to  that  shown  previously  for  spheres 
(Ref.  6)  and  the  comparison  has  been  published  in  summary  form.®*»®2  Comparison  with 
the  sphere  case  reveals  certain  geometrical  differences  in  tiie  backscatteiing  contributions  by 
reverberations.  The  calculation  justified  the  use  of  the  specular  echo  amplitude  to  calibrate 
the  backscatteiing  measurement  snown  in  Fig.  16. 


2.  Focused  Backscatteiing  by  Liquid- HUed  Spnerical  Reflectors. 
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A  previous  ray  analysis  by  Marstun  and  Langley  of  backscattering  by  fluid  spheres 
described  an  amplitude  enhancement  associated  with  axial  focusing  or  glory  scatteiing.^^ 
The  confirmation  by  comparison  with  the  partial  wave  series  (PWS)  had  been  carried  out 
for  spheres  with  an  acoustical  refractive  index  N  =  c/cfluid  <1.  In  the  present  work  the 
analysis  and  comparison  was  extended  by  Marston  and  Sun  to  the  important  practical  case 
of  a  two-chord  glory  ray  for  N  between  V2  and  2.  This  is  the  ray  that  is  relevant  to  focused 
spherical  reflectors  used  by  some  acousticians  for  sonar  calibration  or  for  the  enhancement 
of  target  strength  for  thin  shells  filled  with  a  liquid.  The  comparison  confirms  the  ray 
analysis  for  the  situation  where  shell  effects  are  neglected.®^  The  special  case  of  the  limit 
N  =  2  was  also  analyzed  and  confirmed.  In  that  limit  the  outgoing  wavefront  has  a 
different  shape  and  the  analysis  of  the  focused  scattering  is  different.  Both  cases  are 
summarized  in  the  survey  chapter  (section  4.8  of  publication  B2)  to  illustrate  relationships 
between  wavefront  geometry  and  scattering. 

3.  Diffraction  Catastrophes 

In  research  support  by  Contract  N00014-85-C-0141,  Frederickson  completed  a 
Ph.D.  dissertation  on  transverse  cusp  caustics  produced  by  the  reflection  of  sound  from  a 
curved  metal  surface  in  water.34  Unlike  the  work  on  scattering  by  shells  discussed  in 
Section  A,  the  elastic  response  of  Frederickson's  metal  surface  was  unimportant  to  the 
formation  of  the  caustic:  caustic  formation  was  principally  a  result  of  the  shape  of  the 
reflector.  During  the  final  year  of  present  grant,  aspects  of  Frederickson's  results  were  put 
in  a  form  suitable  for  publication.^^ 

In  research  also  supported  by  Contract  N00014-85-C-0141,  Kaduchak 
photographed  and  studied  the  conditions  for  the  formation  of  a  novel  caustic  in  the 
backscattering  of  light  by  a  penetrable  spheroid.  (The  experiment  was  performed  with  an 
acoustically  levitated  water  drop.)  Preliminary  results  were  sent  to  J.  F.  Nye  of  Bristol 
who  confirmed^^  a  previous  conjecture  by  Marston^  that  the  caustic  classification  was  that 
of  an  Eg  catastrophe.  After  taking  up  a  Ph.D.  project  on  scattering  by  shells, Kaduchak 
returned  to  this  project  briefly  to  examine  his  experimental  results  with  the  aid  of  Nye's 
analysis.15  The  research  serves  to  explore  the  relationship  between  scattering  geometry  and 
caustic  classification  for  a  complicated  situation  not  previously  considered. 

4.  Response  of  Shear  Layers  Separating  from  a  Circular  Cylinder  to  Small- Amplitude 
Rotational  Oscillations. 

J.  Filler  completed  a  Ph.D.  dissertation  on  this  topic  with  support  from  the  previous 
ONR  contract  (N0(X)14-85-C-0141).37  Subsequently  in  1990  while  providing  engineering 


support  to  the  water  tank  facility  used  in  scattering  experiments  (Section  V),  he  put  the 
shear  layer  research  in  publishable  form.J2  experiment  demonstrated  a  novel  fluid- 
structure  interaction  associated  with  flow  past  an  oscillating  cylinder.  The  reader  is  referred 
to  publication  J2  or  to  a  review  of  the  results  by  Griffin.38 

V.  ACOUSTICAL  RADIATION  AND  SCATTERING  FACILITY 

With  the  proposal  for  this  research  came  a  commitment  by  Washington  State 
University  (WSU)  to  provide  laboratory  space  suitable  for  the  water  tank  facility  needed  for 
the  research  described  in  Section  IV  A,  B,  and  C.  The  needed  space  was  remodeled  at  the 
expense  of  WSU  in  1990  while  the  grant  facilitated  the  procurement  of  a  6500  gallon  (12 
ft  diameter  x  8  ft.  deep)  redwood  tank.  A  second  3(X)0  gallon  (8  ft.  diameter  x  8  ft.  deep) 
wooden  tank  was  also  installed  and  the  grant  was  the  principal  source  of  support  for 
positioning  equipment  purchased  for  use  in  these  tanks.  WSU  provided  a  platform  to 
facilitate  access  by  people  and  equipment  to  space  above  the  tanks.  (The  remodeled 
laboratory  space  was  also  used  for  most  of  the  research  on  four-wave  mixing  described  in 
Sec.  IV  C.)  The  large  wooden  tanks  were  essential  to  the  work  and  several  students  and 
other  staff  contributed  to  the  facility  development  and  testing  including  J.  R.  Filler,  D.  H. 
Hughes,  G.  Kaduchak,  T.  J.  Matula,  S.  P.  Parry,  and  H.  J.  Simpson.  The  summary  here 
illustrates  how  the  facility  was  used  in  the  supported  research.  The  facility  has  remained  in 
operation  for  additional  research  on  radiation  and  scattering. 

Figures  36(a)  and  (b)  are  photographs  of  the  facility.  Shown  in  Fig.  36(a)  are  the 
adjacent  wooden  tanks  with  the  larger  of  the  two  tanks  situated  on  the  right  side.  The 
stairway  to  the  platform  above  the  tanks  is  also  visible.  The  tanks  are  supported  by 
pumping  and  water  filtration  systems  that  are  operated  continuously  except  while  data  is 
actually  being  acquired.  Figure  37  illustrates  how  the  large  tank  was  used  to  acquire  the 
backscattering  data  shown  in  Figs.  6-9  where  the  frequencies  used  ranged  from  220  to  430 
kHz.  The  large  tank  was  needed  for  the  experiments  to  prevent  reflections  from  the  walls, 
top,  and  bottom  from  interfering  with  the  signal  of  interest  The  repetition  rate  is  kept 
sufficiently  low  that  reverberations  have  a  negligible  amplitude  at  the  time  of  the  next  ping. 

The  smaller  of  the  wooden  tanks  shown  on  the  left  of  Fig.  36(a)  was  needed  so  that 
the  radiation  and  scattering  experiment  could  be  set  up  and  aligned  separately.  Figure  25 
illustrates  how  the  8  ft  diameter  tank  is  used  in  acoustic  radiation  experiments.  Visible  in 
Fig.  36(b)  is  the  two-dimensional  positioning  system  that  is  supported  above  the  tank,  the 
EMAT  driver  (on  the  left),  and  the  curved  aluminum  plate  that  travels  down  into  the  water. 
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Figure  36  --  (a)  Water  tank  facility  for  radiation  and  scattering  experiments  in  water, 
(b)  Set-up  for  measuring  radiation  from  flexural  wave  tone  bursts  on  plates. 


Figure  37  E) 
backscattering 


riment  in  the  large  tank  for  studying  the  midfirequency  of  the 
ashdl. 
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